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ABSTRACT 


One  of  the  major  problems  with  a  Limited  Area  Model  is  the  introduction  of  error 
from  the  lateral  boundaries.  The  boundary  condition  provides  a  source  of  forcing  to  the 
interior  of  the  model.  This  forcing  typically  originates  from  a  global  model  such  as 
NOGAPS.  The  transition  at  the  boundary  from  one  model  to  another  invariably  produces 
errors.  Traditionally,  the  way  to  minimize  boundary  error  is  to  move  the  boundary  as  far 
away  from  the  area  of  interest  as  possible.  In  this  way,  the  errors  do  not  have  time  to 
infest  the  LAM  with  “bad”  information.  Moving  the  boundary  far  away  from  the  area  of 
interest  increases  the  computational  forecast  load  and  decreases  its  timeliness.  This  study 
looks  at  how  close  the  lateral  boundary  can  be  to  minimize  computational  time  and  still 
maintain  a  forecast  that  is  useful.  It  was  found  that  when  the  entire  inner  COAMPS  nest 
was  analyzed,  the  differences  between  the  control  forecast  and  the  test  forecast  where 
within  the  natural  variability  of  the  control  grid.  It  was  also  found  that  there  where 
localized  areas  within  the  model  domain  that  differed  between  the  control  domain  and  the 
test  domain  by  up  to  20  mb  for  the  sea  level  pressure  after  a  six  day  forecast. 
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I. 


INTRODUCTION 


One  of  the  difficulties  with  Limited  Area  Models  (LAM)  is  the  treatment  of  the 
Lateral  Boundary  Conditions  (LBC).  In  a  global  model,  the  domain  size  encompasses  a 
spherical  globe  and  there  are  no  boundaries.  The  equations  for  the  mathematical  model 
are  closed  in  the  sense  that  there  are  no  lateral  boundary  conditions.  With  such  a  large 
geographical  area  to  cover,  high  resolution  models  cannot  adequately  simulate  mesoscale 
weather  features.  The  model  run  would  be  computationally  intensive  and  operationally 
inhibitive.  LAMs  are  required  to  capture  mesoscale  phenomena  influenced  by  the  local 
terrain  and  local  area  conditions.  One  such  LAM  is  the  Navy’s  Coupled 
Ocean/ Atmosphere  Mesoscale  Prediction  System  (COAMPS).  CO  AMPS  ( Hodur ,  1997), 
obtains  its  LBC  from  the  Navy  Operational  Global  Atmospheric  Prediction  System 
(NOGAPS),  ( Hogan  and  Rosmond,  1991). 

According  to  Warner  et  al.  (1997)  the  LBC  can  influence  the  resultant  solution  of 
a  LAM  based  on  the  treatment  of  the  following  factors.  1)  Interpolation;  the  LBCs  are 
based  on  the  forecast  of  a  coarser  resolution  model.  The  horizontal,  vertical,  and 
temporal  resolution  of  the  boundary  information  is  generally  poorer  than  the  LAM.  This 
necessitates  the  interpolation  of  boundary  information  onto  the  finer  grid  model  domain. 
For  NOGAPS  and  the  COAMPS  81km  nest,  the  grid  sizes  are  nearly  the  same.  Even  so, 
the  geographical  location  of  the  grid  point  may  not  be  coincident,  so  interpolation  would 
still  be  required.  2)  Incorrect  data;  the  model  forecast  provided  from  the  coarser  grid 
domain  may  simply  be  incorrect.  Incorrect  data  could  occur  due  to  the  model  not 
handling  the  true  weather  correctly,  or  because  there  is  a  meteorological  phenomenon  on 
the  boundary  of  the  LAM  that  the  coarse  grid  model  is  not  designed  to  handle.  3)  One¬ 
way  interactions;  longer  wavelengths  generated  or  modified  within  the  LAM  cannot 
feedback  to  the  coarse  grid  model.  In  other  words,  the  solution  to  the  LAM  may  modify 
the  forecast  outcome  of  the  global  model  which,  in  turn,  provides  the  boundary  condition 
of  the  LAM.  With  a  one-way  interaction,  this  feedback  is  not  considered.  4)  Transient 
modes;  specific  LBC  schemes  can  produce  transient  non-meteorological  gravity-inertia 
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modes  on  the  LAM  domain.  These  transient  waves  may  confuse  the  interpretation  of  the 
LAM  forecast.  5)  Parameterizations;  the  coarse  and  fine  grid  models  may  differ  in  their 
handling  of  physical-process  parameterizations.  The  resulting  differences  in  the  solutions 
at  the  boundary  result  in  spurious  gradients  between  the  two  grids. 

Any  errors  generated  on  the  boundary  will  propagate  into  the  interior  of  the  LAM 
along  with  the  mean  flow  of  the  model.  Although  the  magnitude  of  the  influence  of  the 
errors  will  depend  on  factors  such  as  the  strength  of  the  cross-boundary  flow,  strength  of 
the  forcing  at  the  lateral  boundaries  and  the  length  of  the  integration,  errors  in  the  LBC 
will  influence  the  resulting  forecast.  To  mitigate  lateral  boundary  errors,  conventional 
wisdom  has  been  to  move  the  lateral  boundary  as  far  away  from  the  forecast  area  of 
interest  ( Warner  et  al,  1997).  At  this  point  the  question  arises  of  how  close  can  the  outer 
boundaries  be  placed  to  minimize  computational  time,  yet  far  enough  away  to  minimize 
the  adverse  effects  of  the  boundary. 

A.  BACKGROUND 

1.  The  Fundamental  Equations 

The  three  fundamental  laws  that  govern  atmospheric  motion  are  conservation  of 
mass,  conservation  of  momentum  and  conservation  of  energy.  With  the  proper  treatment 
of  these  fundamental  laws,  the  behavior  of  a  parcel  of  air  can  be  described  and  predicted 
to  a  certain  degree.  By  applying  these  laws,  a  relatively  straightforward  set  of 
fundamental  equations  can  be  derived  and  used  to  describe  the  synoptic,  meso-  and 
microscale  fluid  motions  in  the  atmosphere. 

a.  The  Equations  of  Motion 

In  order  to  apply  the  conservation  laws  to  the  atmosphere,  it  is  convenient 
to  start  with  an  analysis  of  the  fundamental  forces  influencing  the  motions  of  air  parcels. 
The  first  of  these  forces  is  the  gravitational  force.  Gravity  is  a  body  force;  the  force 
depends  on  its  mass  and  acts  at  the  center  of  the  body.  Although  gravitational  force 
exists  between  all  bodies,  the  dominant  interaction  for  a  parcel  of  air  is  between  itself  and 
the  earth.  This  interaction  can  be  express  as 
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— —  =  a*  = 

m 


GM( 


\rj 


(0.1) 


where  g*  can  be  considered  a  constant. 


The  second  force  to  consider  is 
the  pressure  gradient  force  (PGF).  The  PGF  is 
an  example  of  a  surface  force  because  it  is 
independent  of  the  mass  and  only  depends  on 
the  surface  characteristics.  The  PGF  arises  due 
to  the  molecular  kinetic  interaction  between  the 
particles  within  a  volume.  To  analyze  the  PGF 

acting  on  a  parcel  of  air  we  consider  an 
infinitesimal  volume  element  SV  =  SxSySz 

centered  at  the  point  xa,  ya,  andz0  seen  in 


(*o .  Vo.zo' 


Figure  1  Pressure  force  acting  on  left 
and  right  face  of  a  unit  cube.  (. From 
Holton  1972 ) 


Figure  1.  If  the  pressure  in  the  center  of  the  volume  is  po,  then  the  pressure  an 

infinitesimal  distance  away  from  the  origin  can  be  expressed  as  a  Taylor  series 
expansion. 


Sp  Sx  5p  8y  Sp  5z  . 

■■  d  H - 1 - 1 - h  higher-order  terms 

0  Sx  2  8y  2  8z  2 


Here  the  non-linear  and  the  higher  order  terms  have  been  ignored.  The  pressure  force 
acting  on  a  volume  of  air  on  side  B  in  the  x  direction  is: 


fBx=- 


(  Sp  Sx 

p  -\ - 

^  °  Sx  2 


SySz 


and  on  side  A  is: 


FAx  “  + 


Sp  Sx 
Sx  2 


SySz 


The  net  pressure  force  acting  on  a  volume  in  the  x  direction  is 
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f,=Fa  +  F.,=~SxSySz 


Sx 


When  dividing  out  the  mass  and  summing  the  pressure  force  in  all  three  orientations  the 
combined  PGF  (per  unit  mass)  becomes 


F  1 

£pgl  =  __LV/? 

m  p 


where  p  is  the  mass  per  unit  volume. 


(0.2) 


The  final  fundamental  force 
acting  on  a  volume  of  air  is  the  frictional 
force.  Like  the  PGF,  frictional  affects  are 
surface  forces  and  do  not  depend  on  the 
mass  of  the  volume.  The  internal  frictional 
force  (viscosity)  arises  due  to  a  fluid’s 

tendency  to  resist  flow.  The  simplepigure  2  Showing  Shearing  Stress  Acting 

example  is  a  fluid  between  two  plates,  Tangential  to  the  Surface  of  a  Unit 

Volume.  ( From  Holton  1972 ) 

one  stationary,  the  other  moving  at  a 

constant  velocity.  The  fluid  adjacent  to  the  stationary  plate  is  at  rest  relative  to  the 
motionless  plate.  The  fluid  adjacent  to  the  moving  plate  is  moving  at  constant  velocity 
with  the  moving  plate.  The  fluid  between  the  plates  must  change  from  a  constant 
velocity  next  to  the  moving  plate  to  no  motion  next  to  the  stationary  plate.  On  the 
molecular  scale,  there  is  a  transfer  of  momentum  of  the  fluid  molecules  between  the 
plates.  By  Newton’s  second  law  there  must  be  a  force  associated  with  this  shearing  stress 
to  maintain  the  upper  plate  in  constant  motion.  The  shearing  stress  is  the  frictional  force 
and  is  proportional  to  the  surface  area  that  the  sheering  stress  is  acting  on. 
Mathematically,  the  sheering  stress  is  defined  as 
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where  zzx  is  the  shearing  stress  in  the  x  direction  owing  to  the  shear  normal  to  the  z 
coordinate,  u  is  the  velocity  in  the  x-direction  and  p  is  a  constant  called  the  dynamic 
viscosity  coefficient. 


In  the  above  example  there  is  no  net  force  acting  on  the  system  of  plates 
and  fluid  since  the  force  acting  to  maintain  the  top  plate  in  motion  is  equal  and  opposite 
to  the  force  required  to  maintain  the  bottom  plate  stationary.  In  the  more  general  case  of 
non-steady  flow  in  compressible  fluid  such  as  air,  a  Taylor  series  expansion  can  be  used 
again  to  determine  an  approximation  to  the  change  of  the  frictional  force  across  a 
differential  volume  of  air.  The  stress  acting  at  a  point  on  the  top  face  of  the  volume  in 
dz  S  z 

Figure  2  is  z2X  + — — —  and  the  stress  acting  at  a  point  on  the  bottom  face  is 

'  dz  2 


drz,  Sz 
dz  2 


Therefore  the  total  stress  acting  on  the  entire  top  and  bottom  face  of 


the  volume  in  Figure  2  is 


+ 


8z' 
dz  2  j 


SySx 


d-^S-fsySx 

dz  2  J 


dz 

—^-SzSySx 

dz 


This  equation  can  be  divided  by  the  mass  of  the  volume  pSxSySz  to 
obtain  the  viscous  force  per  unit  mass  owing  to  vertical  shear  of  the  component  motion  in 
the  x  direction  is 


1  dz  1  d  (  du\ 

- —  = - p — 

p  dz  p  dz\  dz  J 

Here  the  dynamic  viscosity  coefficient  is  nearly  constant  and  can  be 
combined  with  the  density  to  become  the  kinematic  viscosity  coefficient,  v .  The  above 

d2u 

equation  simplifies  to  o — - .  The  resulting  frictional  force  components  per  unit  mass  in 

dz 

the  three  Cartesian  coordinate  directions  are 
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d2u  d2u  d2u 

“  L5x  Qy  dz  . 

„  d2v  d2v  d2v 

„  ["  d2  w  d2  w  d2  w 

Fr=o  — r  + — r  +  — T 
dx2  dy2  dz2 


where  u,  v  and  w  are  the  components  of  velocity  in  the  x,  y  and  z  directions. 
b.  The  Momentum  Equations 

On  the  synoptic  scale  of  motion,  the  reference  frame  for  atmospheric 
motion  is  the  rotating  spherical  surface  of  the  earth.  The  fluid  motion  as  observed  from 
the  rotating  reference  frame  of  the  earth’s  surface  appears  different  than  what  would  be 
observed  from  an  inertial  reference  frame  such  as  the  earth’s  center.  The  change  of  some 
quantity  q  with  respect  to  the  center  of  the  earth  represents  the  absolute  change  in  the 

quantity  and  is  expressed  as  .  Locally,  fixed  on  the  surface  of  the  earth,  that  change 

is  expressed  as  some  local  differentiation,  —  plus  the  advection  of  the  quantity  on  the 

dt 

relative  reference  system,  yVq  .  In  equation  form  this  becomes: 

=  +  V.v  (0.4) 

Dl  dl 

Absolute  change  in  q  Local  change  in  q  Advection  of  q 

In  an  inertial  reference  frame,  Newton’s  second  law  of  motion  states  that 
the  sum  of  the  forces  acting  on  a  body  is  equal  to  the  absolute  change  of  the  momentum 
of  the  body.  In  this  treatment  the  mass  of  the  body  or  parcel  of  air  will  remain  constant 
so  that  any  net  force  acting  on  the  air  parcel  will  induce  a  change  in  its  absolute  velocity, 


^  real  forces  acting  on  parcel  (0.5) 

In  the  inertial  reference  frame  of  the  earth,  there  is  no  apparent  rotational 
affect.  There  is  an  apparent  deviation  in  the  path  of  a  body  or  parcel  when  observed  from 
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the  rotating  reference  frame  relative  to  the  surface  of  the  earth.  Mathematically,  the 
reason  for  this  can  be  shown  by  determining  the  time  rate  of  change  of  a  generic  vector 
quantity  A. 


DA 

Dt 


dA  .  dAv 


dt 


-i  +  - 


dt 


J  +  - 


dA 


dt 


,  rfi  ,  d\  .  dk 
+  A—+AV  —  +  AZ  — 
dt  dt  dt 


dA 

dt 


+  O  x  A 


(0.6) 


where  Q  is  the  rotational  vector  of  the  earth. 

In  order  to  determine  the  time  rate  of  change  of  the  absolute  velocity,  the 
absolute  change  in  the  radial  vector  must  also  be  known. 


Dr 

v„  =  —  =  v  +  t2xr 
Dt 

Dya  d  ,  _  x 

Dt  dt  1 

=  —  +  Qx  v  +  Qx  v  +  fix(Qxr) 

dt  v  ’ 

=—■ +  2Qxv-  QjR 

dt  Coriols  Centrifugal 


Here  R  is  the  position  vector  from  the  axis  of  rotation  toward  the  particle 
on  the  surface  of  the  earth.  It  becomes  clear  that  when  the  reference  frame  shifts  from  the 
absolute  of  the  earth  center  to  the  rotating  surface,  two  apparent  forces  appear;  the 
Coriolis  and  the  Centrifugal.  When  summing  the  forces  on  the  right  hand  side  of 
equation  (0.5),  the  three  real  forces  that  need  to  be  considered  are  the  pressure  gradient 
force,  the  gravitational  force  and  friction  discussed  above.  It  is  convenient  to  combine 
the  gravitational  force  and  the  apparent  centrifugal  force  by  defining  the  gravitational 
force  g  such  that 

g  =  g*+Q2R 

Equation  (0.5)  above  now  becomes 
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Dx  1  1  dr 

- =  -2Q  x  v V p  +  g  4 - 


In  the  familiar  Cartesian 
coordinate  system,  the  velocity  vector  can 
be  represented  by  the  equation 

U  =  ui  +  vj  +  wk  where  the  i ,  j  and  k  are 
orthogonal  unit  vectors.  When  moving 
from  the  inertial  reference  frame  of  the 
earth’s  center  to  a  relative  reference  frame 
of  the  surface,  the  absolute  velocity  vector 
can  be  represented  as 


j  (north) 
s\  i  (east) 
V^k(up) 


d/ d<p  ~  dr  - 

U  =  r  co  s  (p  — —  i  +  r  — —  j  + — k .  F igure  3  Showing  Relationship  between  Earth 

c'  t  tl  at  Centered  Spherical  Coordinates  and  Surface 

Here  i  ,  j  and  k  are  the  new  unit  based  Cartesian  Coordinates. 


vectors  directed  toward  the  east,  the  north  and  up  respectively.  The  distance  from  the 
center  of  the  earth  can  be  represented  by  r  =  a  +  z  where  a  is  the  radius  of  the  earth  and  z 
is  the  height  above  the  surface.  Since  a  »  z ,  r  can  be  replace  by  a  and  is  considered  to 
be  a  constant  when  dealing  with  atmospheric  flow.  For  notational  simplicity,  it  is 
convenient  to  define  x  as  the  distance  in  the  easterly  direction,  y  as  the  distance  in  the 
northerly  direction  and  z  to  point  up  from  the  surface. 

The  velocity  vector  v  of  a  body  on  the  relative  reference  frame  of  the  earth 
is 

v  =  ui  +  vj  +  wk 


Using  equation  (0.6),  the  time  rate  of  change  of  the  velocity  vector  v  is 
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D\  :  du  -  dv  ~  t/w  _ 

- =  1  —  +  j  —  +  k  —  +  Qxv 

Dt  dt  dt  dt 


du  uvtan#?  uw 
dt  r  r 


=  i 

+  j 

+  k 


dv  w2tan<z>  vw 

—  + - —  +  — 

dt  r  r 


dw  u 
dt  r 


2  v2 


Now  equating  the  right  and  left  hand  side  of  equation  (0.7)  and  combining  terms  in  the 
three  coordinate  directions 


du  uv  tan  (p 
dt  r 

dv  u2  tan^ 
dt  r 


uw  _ 
r 


vw 

+  —  = 
r 


1 


dp 


pr  cos  cp  d/ l 


+  2Q( 


vsin^»-  wcos 


<P)  + 


— —  -  2Q  u  sin  (p  +  — - - 

pr  d(p  p  dr 


1  d_lA 

p  dr 


(0.8) 


dw  u2 +v2 
dt  r 


- g  +  2£lu  cos  cp 

p  dr 


To  simplify  the  above  equation  the  following  steps  can  be  taken.  Using 
the  shallow  approximation,  which  states  thatz<sa,  i.e.  the  atmosphere  is  a  thin  layer 
residing  on  the  earth,  and  using  scale  analysis,  the  w  terms  can  be  considered  much 
smaller  than  u  and  v  and  can  be  removed  from  the  first  two  equations  in  (0.8).  Likewise 
the  u  and  v  terms  can  be  removed  from  the  third  equation  in  (0.8).  The  frictional 
contribution  can  be  represented  by  F  and  2Q  sin  cp  is  the  Coriolis  parameter  f.  The 
second  terms  of  the  first  two  equations  can  be  removed  by  assuming  a  local  tangent  plane 
is  being  used.  This  approximation  applies  whenever  the  horizontal  scale  is  much  smaller 
than  the  radius  of  the  earth.  The  final  equations  of  motion  become 
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(0.9) 


Du 

Dt 

Dv 

Dt 

Dw 

Dt 


g 


1  dp 
p  dr 


This  completes  the  derivation  of  the  fundamental  equation  for  the 
conservation  of  momentum  introduced  in  the  beginning  of  the  section.  Equations  (0.9) 
are  referred  to  as  the  equations  of  motion  for  the  atmosphere  on  a  local  plane. 


c.  The  Continuity  Equation 
In  the  next  conservation 


law,  consider  the  mass  contained  in  some 
unit  volume  SV  =  SxSySz .  We  can 
again  represent  the  introduction  of  mass 
on  the  left  side  of  a  cube  as  a  truncated 
Taylor  series  as  shown  in  Fig  4.  The 
amount  of  mass  leaving  the  cube  on  the 
right  side  is  also  represented.  The  change 
in  the  mass  from  some  component  of 
velocity  in  the  x  direction  is 


Figure  4  Showing  the  Transport  of  Mass 
through  a  Unit  Volume.  ( From 


pu  $  yd z  ~  pu  +  -7—  ( pu  )  8  y8z  =  --^-(pu'j  dxSySi 


so  that  the  total  flow  of  mass  from  all  three  directions  becomes 


—  (pu)  +  —  (pv)  +  —  (pw)  8x8yd z 

fix  dy  dz 
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and  the  mass  inflow  per  unit  volume  is  —  V*(/?U).  Since  mass  cannot  be  created  nor 
destroyed  within  the  volume,  any  change  to  the  density  of  the  volume  (— )  would  be  due 


to  the  net  flux  of  mass  of  the  box. 

§ +V.(pU)=0 


(0.10) 


or, 


^  +  /?V.U  +  U»Vp  =  0 


or,  if  the  first  two  terms  are  combined, 


Dp 

Dt 


+  pV»U  =  0 


(0.11) 


Equation  (0.10)  differs  from  equation  (0.1 1)  in  that  the  first  equation  gives 
the  local  rate  of  change  of  the  density  and  the  mass  divergence  and  the  second  equation 
gives  the  absolute  change  in  the  density  of  a  parcel  of  air  following  the  motion  of  the  air 
parcel. 


By  applying  a  synoptic  scale  analysis  to  the  above  continuity  equations,  a 
simplified  form  can  be  used.  Equation  (0. 1 1)  becomes 


dp  dp  dp  dp 

—  +  U  —  +  V—  +  W—+P 

dt  dx  dy  dz 


du  dv  dw 

- 1 - 1 - 

dx  dy  dz 


=  0 


The  density  parameter  p  can  be  approximated  by  the  quantity  p()  +  p' 
where  /?0  is  the  average  density  of  the  volume  and  p  is  the  deviation  from  the  average. 
In  this  way,  time  differentiation  of  the  average  density  term  will  be  zero. 


dp'  dp'  dp'  dp'  dpn  ,  du  dv  dw 

—  +  u  —  +  v—  +  w—  +  w-^-  +  (p  +  p)( —  +  —  +  — )  =  0  (0.12) 

dt  ox  dy  dz  dz  v — v — /  dx  dy  dz 

' - Ts - -  “  '~zZ  ^  z 

~—$P  D  H  D+E  D 

' - ^ '  B  C 

A  F 
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Here,  the  synoptic  scales  used  are 


L 

Horizontal  length  scale 

H 

Vertical  length  scale 

D 

Vertical  length  scale  of  disturbance 

U 

Horizontal  scale  for  advection 

W 

Vertical  scale  for  advection 

In  the  case  of  a  compressible  fluid  such  as  the  atmosphere  the  following 


approximations  can  be  made.  The  density  disturbance  term  Sp  is  much  smaller  than  the 
average  density  of  the  fluid  p() .  Given  that,  term  A  in  equation  (0.12)  is  much  smaller 

than  term  (E)(F)  and  can  be  removed.  Term  B  is  much  smaller  than  term  (D)(G)  and  can 
be  removed.  Finally,  term  E  is  much  smaller  than  term  D  and  can  also  be  removed. 
Equation  (0.12)  is  simplified  to 

w^  +  A>v*v  (0.13) 

dz 

This  equation  will  take  on  different  forms  depending  on  whether  the 
atmospheric  fluid  under  consideration  is  non-hydrostatic  as  in  the  COAMPS  model  (see 
section  3  below)  or  hydrostatic  such  as  the  NOGAPS  model  (see  section  2  below)  where 
the  first  term  in  equation  (0. 13)  can  be  ignored. 

d.  The  Thermodynamic  Equation 

The  law  of  conservation  of  energy  states  that  the  energy  added  to  a  system 
is  equal  to  the  change  in  internal  energy  of  the  system  plus  the  work  done  by  the  system. 
In  the  atmospheric  air  parcel  case  the  internal  energy  is  due  to  the  kinetic  interactions  of 
molecules  within  the  control  volume.  Work  is  done  when  the  boundaries  of  the  control 
volume  is  changed.  The  work  performed  on  the  control  volume  originates  from  the 
external  forces  described  above  and  acting  on  the  surface  area  of  the  control  volume. 
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The  thermodynamic  energy  of  a  control  volume  can  be  represented  in  a 
U*U 

Lagrangian  system  as  p[e-\ — —  ]SV  where  e  represents  the  internal  energy  per  unit 

mass  and  U  is  the  velocity  of  the  control  volume  system.  This  second  term  represents  the 
kinetic  energy  of  the  system.  The  pressure  force  acts  on  the  entire  surface  of  the  control 
volume.  The  rate  at  which  the  pressure  force  does  work  on  the  control  volume  is  the  dot 
product  of  the  force  per  unit  area  and  the  velocity  vector.  Hence,  the  rate  at  which  the 
environment  does  work  on  a  unit  volume  in  the  x  direction  is 


(pu)ASySz ~(pu)BSySz . 


In  this  infinitesimal  control  volume  the  work  performed  on  side  A  is  nearly  equal  to  the 
work  performed  on  side  B.  Therefore 


(pu)B=(pu)A  + 


Sx-\ — 


and  the  net  rate  at  which  pressure  does  work  on  the  control  volume  in  the  x  direction  can 
be  approximated  by 


[(pu)A-(pu)B]SySz 


SV. 


Similar  expressions  can  be  found  for  the  y  and  z  directions.  The  total  rate  at  which 
pressure  does  work  on  the  control  volume  becomes 

-V*(pV)SV. 


As  discussed  earlier,  the  frictional  force  acts  perpendicular  to  the  velocity 
vector  and  therefore  does  not  contribute  to  the  work  performed  on  the  control  volume. 
Likewise,  the  Coriolis  force,  -Q  x  U ,  is  perpendicular  to  the  velocity  and  also  does  not 
contribute  to  the  work.  The  gravity  force  g  can  be  dotted  with  the  velocity  vector  U  when 
directed  downward.  Summing  all  contributions  the  conservation  of  energy  equation, 
when  applied  to  a  control  volume  in  the  Lagrangian  system,  becomes 
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D_ 

~Dt 


(  u.uYI 

p\ 

e  H - 

v  2 

=  -V*(  pV)  8V  +  pg'USV  +  pJ8V 


(0.14) 


Here  J  is  the  rate  of  heating  per  unit  mass  owing  to  radiation,  conduction,  and  latent  heat 
release.  The  last  term  in  equation  (0.14)  represents  the  contribution  to  the  internal  energy 
from  radiation  and  latent  heat  release.  Using  the  chain  rule  (0.14)  becomes 


pSV— 

Dt 


( 


U*U 


e  +  - 
v  2 


f 


+ 


U*U 


e  +  - 
V  2 


D(pSV) 


Dt 


=  -U  •VpSV  -  pV'USV  -  pgwSV  +  pJSV  (0. 1 5) 


where  g  in  equation  (0. 14)  only  acts  parallel  to  the  z  component  of  the  velocity.  From  the 
Lagrangian  form  of  the  continuity  equation  pSV )  =  0  and  the  second  term  on  the  left 


hand  side  of  (0.15)  is  equal  to  zero.  Thus 


De  D 

p - 1 - 

Dt  Dt 


U*U 


=  -U»V/>-/)V»U-yC>gw+  pJ 


(0.16) 


Using  the  Lagrangian  form  of  the  momentum  relation  above  in  equation  (0.7)  without 
friction,  the  following  relation 


P 


can  be  used  to  simplify  (0. 16)  to 


D_ 

~Dt 


1 


\ 


U’U 

v2  j 


=  -U  •Vp-  pgw 


De 

p —  =  -pV»V  +  pJ 
Dt 


(0.17) 


(0.18) 


The  thermal  energy  equation  of  (0.18)  can  be  written  in  more  familiar 
form  by  noting  from  (0. 1 1)  that 

iv.u=-'  Dp -Da 


p 


p  Dt  Dt 


and  that  for  dry  air  the  internal  energy  per  unit  mass  is  given  by  e  =  cDT ,  where  cv  is  the 
specific  heat  at  constant  volume.  Here  a  =  M  p  and  is  called  the  specific  volume.  Now 
the  familiar  form  of  the  thermodynamic  equation  becomes 
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(0.19) 


DT  Da 
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The  above  fundamental  equations  along  with  the  ideal  gas  law  can  be  grouped  together  to 
form  a  complete  set  of  equations  to  describe  the  fluid  motion  of  the  atmosphere. 
Computational  models  can  use  these  equations  to  aid  in  forecasting  and  in  understanding 
how  atmospheric  flows  occur. 


2.  A  Brief  Discussion  on  the  Navy  Operational  Global  Atmospheric 
Prediction  System  (NOGAPS) 

In  the  NOGAPS  model  ( Hogan  and  Rosmond,  1991),  generated  by  Naval 
Research  Laboratories  (NRL)-Monterey,  the  Navy  uses  the  primitive  equations  to 
construct  a  spectral  global  atmospheric  weather  prediction  system.  The  NOGAPS  model 
equations  are  provided  in  Appendix  B.  This  model  is  important  in  providing  guidance 
for  world- wide  naval  operations  and  providing  the  forcing  and  boundary  conditions  for  a 
number  of  atmospheric  and  oceanographic  application  programs.  NOGAPS  provides  the 
real  data  initial  conditions  required  to  set  up  the  first  guess  fields  that  represent  the  initial 
state  of  the  atmosphere  in  the  Coupled  Ocean/ Atmosphere  Mesoscale  Prediction  System 
(COAMPS),  the  mesoscale  atmospheric  model  of  the  United  States  Navy.  In  the 
COAMPS  model  ( Hodur ,  1997),  this  is  referred  to  as  a  ‘NOGAPS  cold  start’  and  is  one 
of  the  two  methods  to  initiate  a  COAMPS  real  data  model  forecast. 

NOGAPS  version  4.0  was  released  in  1998  and  contains  239  spectral  waves  with 
a  0.5-degree  resolution.  There  are  30  sigma  levels  with  approximately  six  levels  below 
850  mb,  depending  on  the  terrain.  NOGAPS  utilizes  a  hybrid  pressure  coordinate 
system.  The  numerical  method  in  NOGAPS  uses  a  horizontal  spectral  differencing, 
second-order  finite  difference  in  the  vertical,  and  central  time  differencing  with  Robert 
semi-  implicit  corrections.  Additional  information  on  NOGAPS  can  be  found  in  Hogan 
and  Rosmond  (1991). 
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3.  The  Basics  of  the  Coupled  Ocean/Atmosphere  Mesoscale  Prediction 
System 

COAMPS  was  implemented  at  the  Fleet  Numerical  Meteorological 
Oceanographic  Center  (FNMOC)  in  October  of  1997.  COAMPS  provides  the  Navy  with 
a  “fmer-scale”  meteorological  representation  of  the  atmosphere  by  decreasing  the  grid 
size  of  the  model  and  introducing  additional  physics  that  can  be  handled  on  a  local  scale 
that  otherwise  could  not  on  the  NOGAPS  global  scale  model.  Numerically,  the  model 
uses  the  primitive  equations  and  includes  nonhydrostatic  effects  on  an  Arakawa  C  grid. 
The  model  equations  used  for  COAMPS  are  provided  in  Appendix  C. 
a.  Grid  Structure 

The  Arakawa  C  grid  was  chosen  due  to  its  better  performance  regarding 
geostrophic  adjustment  and  divergent  flow.  The  grid  consists  of  the  momentum 
components  (u,v,w)  shifted  one-half  grid  interval  from  the  location  of  the  other 
prognostic  variables.  The  grid  is  staggered  in  both  the  horizontal  and  vertical  direction. 
There  are  thirty  sigma  levels  in  the  vertical  in  the  current  operational  configuration.  A 
horizontal  depiction  of  the  Arakawa  C  grid  is  shown  below  for  the  prognostic  variable  n , 
the  dimensionless  pressure. 


Ui~\  ,/+i 

ui,j+ 1 

A+i,/+i 

VMJ 

vu 

v*u 

ui-u 

uiJ 

V  i 

Ui-Uj- 1 

nUi- 1 

A+i,/-i 

b.  COAMPS  Nested  Domains 

The  COAMPS  model  can  be  run  with  up  to  seven  horizontally  nested 
grids.  The  resolution  of  each  inner  consecutive  nest  is  increased  by  a  factor  of  three. 
Every  third  grid  in  the  inner  nest  corresponds  to  a  grid  point  in  the  outer  nest.  The  outer 
grid  is  also  called  the  coarse  nest  for  the  parent/child  nest  pairing.  The  inner  nest  can  also 
be  called  the  fine  nest.  The  current  operational  grid  sizes  used  at  FNMOC  are  8 1 ,  27  and 
9  km  grid  spacing.  Higher  resolution  combinations  can  also  be  used.  FNMOC  will  soon 
be  increasing  the  operational  grid  spacing  to  54,  18,  and  6  km.  The  latter  combination 
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was  the  grid  spacing  used  in  this  study.  COAMPS  provides  for  one-way  and  two-way 
interactions  between  the  coarse  and  fine  grid. 

c.  COAMPS  Integration  Scheme 

COAMPS  uses  second-order  leapfrog  time  differencing  to  integrate  the 
slow  modes  (advective  solution)  using  a  time  step  At  that  depends  on  the  grid  resolution 
and  advective  speeds  to  ensure  computational  stability.  Because  COAMPS  is 
nonhydrostatic,  the  simulations  also  contain  sound  and  gravity  wave  solutions.  Their 
faster  propagation  rates  require  integrations  over  a  smaller  forward  time  step,  and  are 
designated  by  the  number  of  fast  time  steps  per  slow  time  step.  Semi-implicit  vertical 
differencing  and  the  use  of  quasi-compressible  assumptions  allow  COAMPS  to  increase 
the  number  of  fast  time  steps  such  that  the  speed  of  sound  can  be  much  less  than  its 
typical  atmospheric  value  of  -350  m/s. 

d.  Initial  Guess  Field 

The  initial  guess  field  comes  either  from  the  NOGAPS  model  (known  as  a 
‘NOGAPS  cold  start’)  or  from  the  previous  COAMPS  forecast  field.  For  the  NOGAPS 
cold  starts,  global  fields  for  the  pressure  levels  at  0.5  degree  resolution  are  interpolated  to 
the  CO  AMPS  model  grid  points.  The  current  boundary  conditions  are  derived  from  the 
Davies  (1976)  or  Perkey-Kreitzberg  (1976)  treatment  of  NOGAPS  forecast  fields  at  3- 
hour  intervals  and  will  be  discussed  in  greater  detail  below.  Additional  information  on 
COAMPS  can  be  obtained  from  Hodur  (1997). 


B.  CURRENT  TREATMENT  OF  MODEL  LATERAL  BOUNDARY 

CONDITIONS 

1.  The  Lateral  Boundary  Condition  Problem 

Regional  models  require  higher  resolution  grids  and  additional  physical 
parameterizations  in  order  to  better  define  the  atmospheric  flow.  In  turn,  the  larger  grids 
and  parameterization  requires  increased  computational  time  and  inhibits  the  ability  to 
perform  such  a  model  on  the  global  scale  given  operational  time  constraints.  This 
limitation  necessitates  that  regional  models  have  lateral  boundary  around  their  perimeter. 
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The  regional  model  resides  within  a  global  model  such  as  NOGAPS.  The  conditions 
along  the  perimeter  of  the  regional  model  are  derived  from  the  global  model. 

Differences  between  the  two  models  can  lead  to  erroneous  solutions  within  the 
regional  model.  The  higher  resolution  of  the  inner  model  requires  that  the  global  model 
grid  point  data  be  interpolated  to  provide  data  to  the  unmatched  inner  grid  points  along 
the  boundary.  In  effect,  the  interpolation  is  generating  mesoscale  data  where  previously 
there  was  none.  This  could  lead  to  erroneous  data  propagating  within  the  regional  model. 
Another  problem  is  the  possible  inconsistent  handling  of  physical  parametric  processes 
between  the  global  and  region  models.  Different  methods  of  parameterization,  such  as 
cumulus  convection,  can  lead  to  erroneous  gradients  in  the  boundary  region  in  which 
error  will  propagate  (Warner  et  al.  1997).  Temporal  resolution  can  also  be  a  problem. 
COAMPS  uses  the  12  hour  NOGAPS  forecast  to  generate  the  COAMPS  coarse  nest 
lateral  boundary  conditions.  Although  the  12  hour  forecast  is  generally  considered  quite 
good  (FNMOC,  2004),  small  discrepancies  can  cause  errors  to  propagate  through  the 
model. 

Proper  handling  of  the  lateral  boundary  conditions  is  required  to  limit  the 
propagation  of  errors  through  regional  models.  The  most  recent  NOGAPS  forecast  is 
used  to  generate  the  COAMPS  outer  nest  lateral  boundary  condition  for  a  real  data 
simulation.  The  NOGAPS  fields  are  horizontally  interpolated  using  a  bicubic  spline  to 
the  COAMPS  coarse  grid,  then  interpolated  in  the  vertical  to  the  model  vertical 
coordinates.  COAMPS  can  either  use  the  Perkey-Kreitzberg  (1976),  which  specifies  the 
boundaries  according  to  the  NOGAPS  tendencies,  or  the  Davies  (1976)  method,  which 
uses  the  actual  NOGAPS  fields  to  determine  the  lateral  boundary  conditions. 

a.  Perkey-Kreitzberg  Method 

In  the  Perkey-Kreitzberg  method,  the  boundary  conditions  consist  of 
large-scale  time- varying  tendencies  provided  by  the  global  model  linearly  combined  with 
model-calculated  tendencies.  The  prediction  of  any  dependent  variable  %  can  be  written 
as: 
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where  the  subscripts  n  and  p  denote  the  “new”  value  after  the  boundary  condition  and  the 
“previous”  value  at  a  former  time.  The  Is  denotes  the  large-scale  specified  tendency  of  x, 
and  m,  the  model-calculated  tendency.  The  values  of  the  weighting  coefficients  W  are 
constructed  such  that  the  interior  grid  points  are  of  value  unity,  the  boundary  grid  points 
are  zero,  and  there  is  a  transition  “buffer”  layer  between  the  two.  In  COAMPS,  the 
transition  is  seven  grid  points  wide.  Thus,  the  outer  grid  points  of  the  LAM  are 
completely  determined  by  the  global  model,  the  interior  points  are  derived  from  the 
regional  model,  and  there  is  a  limited  transition  zone  between  the  two. 

The  Perkey-Kreitzberg  model  has  a  tendency  to  transfer  the  energy  of 
interior  generated  long  waves  to  the  short  wave  region  (Perkey  and  Kreitzberg,  1976). 
The  energetic  short  waves  are  reduced  by  introducing  a  period  filter  in  the  boundary 
region.  In  their  paper,  Perkey  and  Kreitzberg  demonstrate  that  an  exterior  generated 
wave  will  propagate  through  a  regional  model  with  little  deformation.  Interior  waves 
generated  by  physical  processes  within  the  model  will  dissipate  within  the  transition 
zone,  effectively  “leaving”  the  model  region  with  little  reflection  at  the  boundary.  This 
method  proved  to  be  an  effective  way  to  handle  the  boundary  problem  between  global 
and  regional  models. 

b.  Davies  Method 

The  Davies  approach  is  to  construct  a  boundary  formulation  that  is 
consonant  with  the  governing  equations  but  also  alleviates  the  effect  of  the  inaccurate 
boundary  data.  In  this  way,  alternate  governing  equations  are  used  that  reduce  the 
sensitivity  to  boundary  conditions  and  reduce  the  generation  of  noise  on  the  boundary. 
For  example,  in  the  case  of  a  linear,  inviscid,  compressible  flow  confined  between  two 
boundaries  (flow  is  independent  in  the  y-direction),  the  governing  equations  can  be 
written  as: 

ut+Uux  +  wzU  -  fv  =  -px-K(u-u)  (0-21) 
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vt+Uvx  +  f(u  +  U)  =  E-K(v-v)  +  F[v-v] 


(0.22) 


Pi  +  Upx  +  c2ux  =  -c2(wz  +  wS)-K(p-  p).  (0.23) 


Here,  E  ( E=fU)  is  the  external  forcing,  and  the  function  K  and  the  functional  F  contain 
the  mathematical  treatment  of  the  boundary  conditions.  The  mean  flow  is  U  and  an 
estimate  of  the  time  development  of  the  flow  field  in  the  neighborhood  of  the  boundary  is 
denoted  by  the  over  barred  fields  u,  v,  p .  The  specified  function  K  =  K(x)  >  0  is 
continuous  and  non-zero  only  in  the  vicinity  of  the  boundary.  The  functional  F  is  defined 

by  F\v -  v]  =  J  Kx(v -  v)dx  .  In  effect,  K  and  F  act  to  reduce  the  prognostic  variables 

near  the  boundary  that  are  different  than  the  prescribed  values  generated  from  the  global 
model.  The  relaxation  coefficients  smooth  the  fields  in  the  location  of  the  boundary. 

For  the  Davies  methodology,  the  K  term  acts  as  the  weighting  functions 
much  like  the  W  term  discussed  above  for  the  Perkey-Kreitzberg.  The  difference  is  how 
the  weighting  factors  are  used.  For  the  Perkey-Kreitzberg,  as  can  be  seen  in  equation 
1.20,  the  new  prediction  for  the  dependent  variable  within  the  boundary  zone  is 
determined  by  a  time  integrated  transition  of  the  outer  model  to  the  inner  model.  The 
Davies  model  uses  the  boundary  zone  weighting  functions  within  the  governing 
equations,  which  include  both  time  and  spatial  dependent  variables.  The  Davies 
technique  implements  the  boundary  correction  while  solving  the  governing  equations. 
The  Perkey-Kreitzberg  technique  implements  the  boundary  correction  after  solving  each 
time  step  of  the  governing  equations. 

Although  either  the  Perkey-Kreitzberg  or  the  Davies  boundary  treatment 
can  be  selected  in  COAMPS,  the  Davies  method  is  the  only  one  that  has  been  rigorously 
tested  by  NRL  -  Monterey,  and  used  operationally  at  FNMOC  (Doyle,  2004). 
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II.  THE  EXPERIMENT 


A.  EXPERIMENT  SET  UP 

To  analyze  the  effect  of  the  lateral  boundary  on  the  inner  nest  of  a  CO  AMPS 
forecast,  the  following  experiment  was  constructed.  Four  sets  of  COAMPS  forecast  data 
runs  were  performed.  The  coarse  and  fine  domain  of  the  control  run  corresponded  to  the 
first  and  second  nest  of  the  COAMPS  European  Operational  Area  run  at  FNMOC.  The 
coarse  grid  for  the  remaining  three  test  runs  expanded  geographically  to  the  west  and  the 
north  of  the  control  run  coarse  grid.  The  fine  grid  for  the  three  test  runs  had  roughly  the 
same  geographical  position  as  the  control  run  fine  grid.  The  latitude  and  longitude  grid 
points  did  not  match  up  exactly  for  all  runs.  This  slight  difference  in  grid  point 
positioning  will  be  discussed  in  more  detail  below.  The  inner  and  outer  domain 
boundaries  can  be  seen  in  Figure  5. 

Research  on  the  COAMPS  lateral  boundaries  was  performed  in  conjunction  with 
efforts  at  FNMOC  to  study  methods  to  improve  operational  handling  of  the  lateral 
boundaries.  Data  for  this  analysis  was  acquired  from  FNMOCs  project,  “COAMPS 
Operational  Configuration  Testing”.  A  copy  of  the  FNMOC  proposal  is  provided  in 
Appendix  D.  FNMOC  research  was  conducted  over  a  three  month  period  from  15 
December  2003  until  15  March  2004.  COAMPS  forecasts  carried  out  by  FNMOC 
extended  to  48  hours.  For  this  study,  a  request  was  granted  to  extend  the  forecast  fields 
out  to  144  hours  by  rerunning  the  forecasts  using  the  CO  AMPS  restart  data  fields  for  the 
selected  model  dates  of  interest. 

Two  forecast  periods  of  interest  were  chosen.  The  study  set  consisted  of  the 
COAMPS  Analysis  and  the  following  144  hours  of  forecast  fields  at  six  hour  intervals. 
Time  periods  containing  an  “active”  and  an  “inactive”  period  were  chosen  for  study. 
Time  periods  were  chosen  based  on  the  strength  of  the  300  mb  jet  running  zonally 
through  the  upwind  portion  of  the  FNMOC  European  Operational  Area.  The  active  time 
period  contained  a  strong  jet  located  on  the  upwind  side  of  the  European  Operational 
Area  for  a  six  day  period  based  on  twelve  hour  NOGAPS  Analysis  fields  (see  Figure  6  in 
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Appendix  A).  The  model  ran  chosen  for  the  active  period  was  07  January  2004,  0000 
UTC.  For  the  inactive  period,  no  sustained  zonal  jet  was  present  within  the  FNMOC 
European  Operational  Area  (see  Figure  7  in  Appendix  A).  The  model  ran  chosen  for  the 
inactive  period  was  04  March  2004,  0000  UTC. 

For  the  time  periods  stated,  four  sets  of  CO  AMPS  model  runs  were  performed. 
Each  ran  consisted  of  an  inner  fine  nest  at  18  km  resolution  with  the  same  given 
geographical  location.  Each  had  a  coarse  outer  grid  nest  at  54  km  resolution,  with  an 
expanding  geographical 

location.  The  control  data  ran 
is  identical  with  the  inner  and 
outer  grids  used  by  FNMOC 
for  the  operational  European 
grid.  The  control  model  ran 
will  also  be  identified  as 
Runl.  The  remaining  model 
runs  were  identified  as  Run3a, 

Run3b  and  Run3c.  Run3c  had 
the  largest  coarse  grid  and 

Figure  5  Grid  Locations  for  the  Inner  and  Outer 
Run3a  had  the  coarse  grid  Grids  for  Runl,  Run3a,  Run3b  and  Run3c. 

closest  to  the  control  outer  Inner  §rid  for  a11  four  runs  shown  in  Sreen- 

Outer  grid  for  control  ran  shown  in  red.  Outer 
grid.  The  geographical  grid  for  Run3a/b/c  shown  as  dotted,  dot-dashed 

locations  for  each  outer  grid  and  dashed  respectively. 

are  shown  in  Figure  5.  A  spread  sheet  containing  model  set  up  information  derived  from 
the  COAMPS  datahd  fields  are  provided  in  Appendix  E.  The  datahd  file  is  the  file 
generated  from  the  CO  AMPS  Analysis  program  that  indicates  the  domain  configuration 
and  program  settings  used  in  the  particular  ran. 

1.  Data  Field  Mismatch 

The  boundary  locations  for  the  coarse  grids  shown  in  Figure  5  were  determined 
from  the  outer  grid  points  of  the  latitude  and  longitude  fields  from  the  CO  AMPS  runs. 
The  values  contained  in  the  lat/long  fields  correspond  to  the  locations  of  the  data 
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contained  in  the  COAMPS  data  fields.  As  can  be  seen  in  Figure  5,  information  for  Run3c 
has  been  shifted  toward  the  northwest  compared  to  the  other  three  data  runs.  As  a 
consequence,  the  inner  grid  for  Run3c  has  also  been  shifted  (not  shown  in  Figure  5).  The 
mismatch  was  caused  by  the  improper  setup  of  the  Run3c  domain. 

To  get  the  most  accurate  comparison  of  data  points,  the  grid  points  used  in  the 
study  were  compared  based  on  their  geographical  locations,  and  not  necessarily  grid 
point  number.  The  only  shift  required  was  for  Run3c.  In  effect,  the  top  six  rows  of 
horizontal  data  points  of  Run3c  and  the  bottom  six  rows  of  horizontal  data  points  of 
Runl/3a/3b  were  not  included  for  comparison.  As  listed  in  Appendix  E,  this  reduced  the 
number  of  north/south  grid  points  from  196  to  191.  The  total  reduction  in  data  points 
from  57,820  to  56,345  reduces  the  amount  of  analyzed  data  by  2.55%. 

After  the  readjustment  of  grid  points  for  data  field  comparisons,  the  maximum 
deviation  between  the  control  run  latitude/longitude  grid  points  and  the  latitude/longitude 
points  associated  with  the  data  runs  for  Run3a/b/c  was  determined  to  be  small  enough  to 
allow  for  sufficient  comparison.  The  maximum  deviations  in  degrees  and  in  kilometers 
are  listed  in  Table  1  below. 


Run3a 

Run3b 

Run3c 

Max  lat  deviation 

0.0020° 

0.22  km 

0.0082° 

0.90  km 

0.0032° 

0.35  km 

Max  long  deviation 

0.0057° 

0.63  km 

0.0121° 

1.33  km 

0.0059° 

0.65  km 

Table  1  Max  Lat/Lon  Deviation 


2.  Comparison  Method 

To  get  a  visual  picture  of  how  the  separation  distance  between  the  coarse  grid  and 
the  fine  grid  lateral  boundaries  change  the  forecast  outcome,  contours  of  each  data  run 
were  plotted  on  the  same  figure.  The  contour  lines  had  the  following  coloring 
convention;  red  for  the  control  run,  blue  for  Run3a,  green  for  Run3b  and  magenta  for 
Run3c.  To  determine  whether  differences  in  the  contour  lines  were  significant,  the 
following  method  was  used.  The  standard  deviation  for  each  forecast  hour  of  the  control 
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run  was  calculated  using  MATLAB.  Then  MATLAB  was  used  to  determine  how  the  test 
runs  differed  from  the  control  run.  The  control  run  field  was  subtracted  from  the  test  run 
field  point  for  point  over  the  entire  inner  nest  domain.  Then  the  root  mean  square  (RMS) 
difference  of  the  comparison  field  was  calculated.  If  the  RMS  difference  of  the 
comparison  field  was  below  the  standard  deviation  of  the  control  field  alone,  then  it  was 
assumed  that  the  differences  observed  in  the  contour  fields  were  within  the  level  of 
natural  variation  of  the  control  field,  and  therefore  insignificant.  As  will  be  shown,  it 
proved  sufficient  to  compare  the  RMS  differences  with  the  minimum  standard  deviation 
of  the  control  field. 

Runl,  with  the  smallest  outer  domain,  was  used  as  the  control  instead  of  Run3c 
with  the  largest  outer  domain,  as  typically  would  have  been  done.  It  is  assumed  that  the 
run  with  the  largest  outer  boundary  will  have  the  least  adverse  effects  from  the  boundary 
and  therefore  be  the  closest  to  the  true  solution  (Warner  et  al,  1997).  Since  this  work  was 
performed  in  conjunction  with  FNMOC  research,  and  in  their  research  Runl  was  used  as 
the  control  (see  Appendix  D),  the  same  convention  was  used  here.  As  a  consequence, 
computed  differences  between  the  control  (Runl)  and  the  test  runs  (Run3a/b/c)  will 
appear  to  increase  more  quickly  for  the  largest  outermost  domain,  and  more  slowly  for 
the  smallest  outermost  domain. 

A  smaller  outer  domain  would  mean  the  outer  and  inner  boundaries  are  closer 
together.  This  would  imply  that  adverse  boundary  effect  would  populate  the  inner  domain 
more  quickly  and  grow  larger  than  for  a  model  run  where  the  outer  domain  was  larger. 
Here,  the  comparison  is  performed  with  the  model  run  with  the  smallest  outer  domain,  or 
the  one  least  likely  to  be  correct  when  compared  to  truth.  Since  Run3c  has  the  largest 
outer  domain  and  is  most  likely  to  be  correct  when  compared  to  truth,  its  RMS 
differences  grow  the  fastest  when  compared  to  Runl,  the  control. 
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B.  ANALYSIS  OF  RESULTS 

1.  Active  Period 

Figure  8  in  Appendix  A  shows  the  sea  level  pressure  contour  lines  for  the  active 
periods.  The  figure  includes  the  COAMPS  forecasts  in  24  hour  increments  out  to  144 
hours.  The  24  hour  forecast  sea  level  pressure  fields  were  missing  from  the  FNMOC  data 
set,  and  therefore  were  not  included  in  the  analysis.  The  contour  lines  depicted  are  the 
980,  1000  and  1020  mb  contours. 

For  the  SLP  COAMPS  Analysis  (‘000’  forecast  hour)  in  the  upper  left  comer  of 
Figure  8,  the  contour  lines  appear  to  nearly  coincide,  although  there  are  slight  differences. 
The  1020  mb  contour  located  off  the  coast  of  Portugal  splits  into  two  groups.  Contour 
lines  for  Runl  and  Run3a  follow  a  Northeast-Southwest  direction,  and  then  cut  in  to  an 
East-West  direction  before  rejoining  the  other  two  test  run  contours.  The  same  contours 
for  test  runs  3b  and  3c  maintain  a  North-South  orientation  in  the  same  general  location. 
Since  all  runs  had  the  same  settings  and  input  parameters,  slight  variations  in  the  contour 
lines  such  as  these  are  solely  due  to  the  geographical  placement  of  the  outer  grid.  Slight 
variations  in  the  COAMPS  Analysis  occur  solely  due  to  the  geographic  configuration  of 
the  grid  domain. 

In  Figure  9,  a  horizontal  depiction  of  the  difference  between  the  control  and 
Run3c  for  the  COAMPS  Analysis,  ‘000’  forecast  for  the  active  time  period  can  be  seen. 
In  Figure  9,  blue  areas  represent  areas  where  the  test  ran  sea  level  pressure  is  below  the 
control  ran  sea  level  pressure.  Yellow/red  areas  represent  areas  where  the  test  ran 
pressure  is  greater  than  the  control  ran.  Control  ran  contour  lines  have  been  added  for 
reference.  As  can  be  seen  by  the  colorbar  in  Figure  9,  the  difference  between  test  Run3c 
and  the  control  ran  is  slight.  There  are  differences  between  the  control  and  the  test  ran 
COAMPS  Analysis  that  are  solely  due  to  the  geographical  domain  set  up,  but  these 
differences  are  slight,  and  have  little  effect  on  the  sea  level  pressure  contour  lines. 

As  the  COAMPS  forecast  progresses,  the  variability  in  the  SLP  contour  lines 
increase  as  can  be  seen  in  Figures  8.  For  the  active  period  the  variability  is  greatest  in  the 
northwest  region  at  the  120  forecast  hour.  Again,  looking  at  the  horizontal  difference 
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depictions  between  the  control  run  and  Run3c  for  the  active  period,  variations  can  be 
seen.  The  120  and  144  SLP  fields  are  shown  in  Figures  10  and  11  respectively. 
Variations  of  up  to  20  mb  differences  can  be  seen  between  the  control  run  and  Run3c. 

In  Figure  12,  the  standard  deviation  of  the  control  run  is  plotted  over  the  forecast 
period  for  both  the  active  and  the  inactive  period.  The  minimum  standard  deviation  for 
the  active  period  occurred  just  prior  to  the  72  hour  forecast  and  reached  10.6  mb.  The 
minimum  standard  deviation  for  the  inactive  period  occurred  at  the  90  hour  forecast  and 
reached  6. 1  mb. 

In  Figures  13,  the  minimum  standard  deviations  were  compared  to  the  RMS 
differences  between  the  control  run  and  the  individual  test  runs.  The  maximum  RMS 
differences  occurring  at  the  end  of  the  forecast  period  are  far  below  that  of  the  standard 
deviation  of  the  control  run.  For  the  active  time  period  in  Figure  13,  the  maximum  RMS 
was  about  5.5  mb,  which  is  approximately  half  the  value  of  the  minimum  control  run 
standard  deviation.  So  when  the  entire  grid  domain  is  under  consideration,  the 
differences  observed  between  the  test  runs  and  the  control  run  are  within  the  natural 
variability  of  the  control  run  alone. 

As  can  be  seen  in  Figures  10  and  11,  there  are  isolated  areas  where  the  SLP 
difference  between  the  control  and  Run3c  are  as  large  as  20  mb.  For  forecasting 
purposes  a  20  mb  error  could  mean  the  difference  of  predicting  an  active  cyclone  or  not. 
For  instance,  in  Figure  1 1  the  control  run  predicts  a  980  mb  pressure  low  over  central 
Europe.  In  the  same  location,  Run3c  predicts  a  pressure  of  15  to  20  mb  greater, 
equivalent  to  the  surrounding  area  and  therefore  no  low  pressure  system  was  predicted. 
Similarly,  on  the  northwest  edge  of  the  domain,  Run3c  predicts  the  980  mb  contour  to  be 
much  more  extensive  than  the  control.  The  control  run  closes  the  980  mb  contour  with 
the  edge  of  the  domain  to  the  northeast  of  the  British  Island  whereas  Run3c  extends  the 
980  mb  contour  farther  west. 

Even  though  differences  between  the  test  runs  and  the  control  were  found  to  be 
minimal  when  the  entire  domain  was  considered,  there  were  isolated  areas  where 
significant  differences  occur.  From  a  forecaster’s  perspective,  the  “interesting”  weather 
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observed  at  the  surface  is  forced  by  the  presence  of  a  low  pressure  system.  As  the 
cyclonic  motion  wraps  around  and  into  the  low  pressure  system,  high  winds  develop  at 
the  surface.  Higher  winds  drive  the  surface  wave  heights.  For  a  naval  vessel,  high  seas 
and  winds  can  cause,  at  best,  uncomfortable  living  conditions  onboard  ship  and,  at  worst, 
disastrous  damage  to  ship  and  crew.  To  avoid  adverse  weather  conditions,  naval  vessels 
rely  on  weather  forecasters  to  predict  the  location  and  strength  of  low  pressure  systems. 
At  times,  operational  restraints  require  naval  vessels  to  traverse  in  close  proximity  to  low 
pressure  systems,  so  accurate  and  timely  forecasts  are  a  must.  Operational  commanders 
require  the  most  accurate  information,  including  weather  information,  to  plan  effectively. 
Differences  in  forecasted  cyclonic  positioning  and  strength  due  to  differences  in  model 
setups  could  change  the  operational  effectiveness  of  the  naval  vessel. 

To  see  where  the  domain  size  begins  to  affect  the  predicted  outcome  of  the 
forecast,  a  subset  of  the  domain  grid  points  was  used.  The  forecaster  is  concerned  with 
low  pressure  systems,  therefore  for  the  active  case,  the  subset  of  domain  points  below  the 
990  mb  pressure  level  was  considered  for  reevaluation.  The  same  method  of  comparing 
the  STD  of  the  control  with  the  RMS  of  each  run  was  used  to  reevaluate  the  990  mb 
subset.  Here,  the  total  number  of  points  changed  with  each  forecast  time  step,  since  the 
990  mb  contour  changes  with  each  time  step.  Also,  all  four  data  runs  were  used  to 
determine  the  subset  of  data  points  used  to  calculate  the  STD  and  RMS  values  at  each 
time  step. 

Figure  14  shows  the  STD  and  RMS  values  from  the  active  period  sea  level 
pressure  below  990  mb.  The  STD  of  the  control  run  is  in  red.  The  RMS  difference 
between  the  control  run  and  Run3a/b/c  are  given  in  blue,  green  and  magenta, 
respectively.  Below  Figure  14  is  a  table  stating  the  percentage  of  domain  points  having  a 
sea  level  pressure  below  990  mb.  The  percentage  was  determined  after  the  correction 
was  made  for  the  domain  mismatch  discussed  in  the  previous  section.  By  the  72  hour 
forecast,  the  control  run  STD  fell  to  half  its  initial  value  and  the  RMS  values  began  to 
increase.  At  the  78  hour  forecast,  Run3c  RMS  values  exceed  the  control  run  STD. 
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The  Petterssen  Development  Equation  can  be  used  as  a  model  of  the  development 
of  synoptic  scale  extra-tropical  low  pressure  systems.  It  can  be  expressed  as 


v  dt  j  p„ 
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where  the  left  hand  side  represents  the  change  in  the  geostrophic  vorticity  at  the  surface. 


If 


v  dt  jj. 


>0,  cyclogenesis  is  occurring  and  a  low  pressure  system  initiates  and/or 


strengthens.  If 


v  dt  j p 


<0,  cyclolysis  is  occurring  and  a  low  pressure  system  fills  or  is 


inhibited.  The  first  term  on  the  right  hand  side  represents  the  advection  of  vorticity  at  the 
level  of  non-divergence.  Positive  vorticity  advection  (PVA)  at  the  level  of  nondivergense 
contributes  to  cyclogenesis  by  promoting  cyclonic  motion  within  the  air  column  and 
inducing  vertical  motion.  The  second  term  on  the  right  hand  side  represents  the  thermal 
advection.  By  making  simplifying  assumptions  the  warm  air  advection  can  be 
represented  as  a  sine  wave.  After  the  assumptions  have  been  made,  the  Petterssen 
Development  Equation  predicts  that  warm  air  advection  contributes  to  cyclogenesis.  The 
third  term  is  the  braking  term  which  inhibits  development  by  cooling  the  column  due  to 
adiabatic  cooling  of  rising  air.  The  last  expression  on  the  right  hand  side  is  the  diabatic 
heating  term.  The  two  major  development  terms  are  the  vorticity  advection  and  thermal 
advection.  A  closer  study  of  these  two  terms  was  performed  to  get  a  sense  of  how  the 
COAMPS  outer  boundary  impacts  cyclonic  development  using  the  Petterssen 
Development  Equation. 

a.  Case  1 

The  84  hour  forecast  of  the  thermal  advection  is  shown  in  Figure  15  for 
the  active  period.  The  SLP  for  the  equivalent  time  period  is  shown  in  Figure  16,  as  well 
as  the  6  hour  prior  forecast.  As  can  be  seen,  there  is  a  large  area  of  positive  thermal 
advection  over  Central  Europe.  The  fact  that  this  advection  is  stronger  for  the  control 
run,  seen  as  a  deeper  shape  of  yellow,  may  be  contributing  to  the  forecasted  small  low 
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pressure  system  for  the  control  run,  seen  as  the  1000  mb  contour  over  the  Czech 
Republic,  that  is  not  present  for  Run3c,  seen  as  a  yellow  area  in  Figure  16  over  the  same 
location. 

b.  Case  2 

In  the  96  hour  forecast,  large  areas  of  vorticity  advection  for  the  control 
run  active  period  can  be  seen  in  Figure  17.  In  Figure  17,  areas  of  positive  vorticity 
advection  (PVA)  are  shaded  yellow/red.  No  substantial  PVA  or  negative  vorticity 
advection  (NVA)  was  present  in  the  6  hour  forecasts  prior  to  and  following  the  96  hour 
forecast  (figures  not  shown.)  One  area  of  PVA  in  the  96  hour  forecast  is  associated  with 
the  500  mb  trough  extending  over  the  Mediterranean  Sea.  The  influence  of  the  PVA  can 
be  seen  in  the  sea  level  pressure  forecasts  in  Figures  18  and  19.  Figures  18  and  19  shows 
the  progression  of  the  SLP  from  the  90  hour  forecast  to  the  108  hour  forecast.  At  the  90 
hour  frame,  a  previously  developed  relative  low  pressure  system  resides  of  the  country  of 
Greece.  By  the  96  hour  forecast,  the  low  pressure  system  deepened  in  the  control  run 
seen  as  an  expansion  of  the  1000  mb  contour  line.  For  the  next  12  hours,  after  the  PVA 
has  expired,  there  was  no  further  deepening  of  the  low  pressure  system.  In  contrast,  the 
low  pressure  system  for  Run3c  begins  to  fill  between  the  90  and  96  hour  forecast.  This 
can  be  seen  as  a  progression  from  blue  shading,  indicating  the  test  field  below  the  control 
field,  to  yellow  shading,  indicating  the  test  field  above  the  control  field.  By  the  108  hour 
forecast  it  appears  the  control  run  low  pressure  system  extends  from  the  Aegean  Sea  up 
into  Romania.  The  low  pressure  system  for  Run3c  has  completely  filled  over  the  Aegean 
Sea,  and  has  moved  exclusively  over  the  European  Continent  over  Ukraine. 

2.  Inactive  Period 

The  comparison  all  four  COAMPS  Analysis  (‘000’  forecast)  runs  is  show  in 
Figure  20  for  the  inactive  period.  Like  the  active  period  seen  in  Figure  8,  the  contour 
lines  are  closely  packed,  only  indicating  slight  differences  between  the  four  runs.  The 
largest  difference  seen  in  Figure  20  is  the  cut  off  1020  mb  low  over  the  Atlantic.  Figure 
21  shows  the  horizontal  depiction  of  the  difference  between  the  control  run  with  the 
smallest  outer  domain  and  Run3c  with  the  largest  outer  domain.  Here  the  two  fields 
agree  to  within  2  mb  over  the  entire  inner  nest  domain. 
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As  the  forecast  progresses,  the  contour  lines  begin  to  separate  as  can  be  seen  in 
Figure  20.  Like  the  active  period,  the  largest  contour  spreads  seem  to  occur  over  the 
Atlantic  at  the  upwind  side  of  the  inner  nest  domain.  There  is  also  a  large  disagreement 
by  the  120  hour  forecast  on  the  position  and  extent  of  the  1020  mb  relative  low  pressure 
system  in  center  of  the  model  domain.  The  smaller  outer  grid  domains  of  Run3a  and  the 
control  maintain  the  1020  mb  system  over  the  Mediterranean  Sea  whereas  the  larger 
outer  grid  domains  of  Run3b  and  Run3c  extend  the  1020  mb  contour  up  into  lower 
Europe.  Figure  22  shows  the  horizontal  depiction  of  the  numerical  difference  between 
the  sea  level  pressures  of  the  two  extreme  cases  of  outer  domain  size.  The  darkened  blue 
area  over  Ukraine  signifies  that  Run3c  SLP  was  below  the  control  run  in  this  area.  The 
yellow/red  area  over  the  1020  mb  control  run  contour  indicates  that  the  SLP  was  much 
higher  in  this  region  for  Run3c.  The  low  pressure  system  for  Run3c  was  shifted  further 
to  the  northeast  of  the  control  run  low  pressure  system. 

As  with  the  active  period,  a  comparison  was  made  of  the  RMS  difference  between 
the  test  runs  and  the  control  run,  and  the  standard  deviation  of  the  control  run  alone. 
Figure  23  shows  that  the  RMS  differences  during  the  inactive  period  were  below  the 
minimum  value  of  the  control  run  standard  deviation.  The  minimum  standard  deviation 
of  the  control  run  during  the  inactive  period  was  smaller  than  the  minimum  standard 
deviation  of  the  control  run  during  the  active  period.  Therefore  the  maximum  RMS 
differences  were  closer  to  the  minimum  standard  deviation  of  the  control  run.  Still,  this 
shows  that  when  the  entire  inner  domain  field  is  considered,  the  differences  observed  in 
the  test  runs  are  within  the  natural  variability  of  the  control  run,  and  are  therefore 
insignificant. 

Like  the  active  period,  there  are  locations  within  the  domain  were  deviations  are 
large  and  exceed  the  level  of  natural  variability  of  the  control  run.  Again,  taking  the 
forecasters  perspective,  forecasted  differences  of  low  pressure  systems  can  drastically 
change  operational  priorities  for  the  U.S.  Navy.  Using  the  same  analysis  discussed  in  the 
active  section,  the  SLP  contour  below  which  RMS  differences  exceed  the  control  run 
standard  deviation  was  found  to  be  1010  mb,  as  seen  in  Figure  24.  In  the  1010  mb  subset 
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of  domain  grid  points  shown  in  Figure  24,  the  standard  deviation  of  the  control  run  was 
only  3  mb  and  the  RMS  differences  of  the  test  runs  exceed  this  value  around  96  hours. 
By  144  hours,  the  RMS  differences  are  approximately  the  same  value  as  the  control 
maximum  standard  deviation  value.  At  this  point  the  differences  between  the  model  runs 
become  significant. 

The  Petterssen  Development  Equation  was  again  used  to  see  how  differences  in 
the  upper  atmosphere  changed  the  level  of  surface  cyclonic  development. 

a.  Case  3 

Figure  25  shows  the  advection  of  relative  vorticity  for  the  inactive  period 
at  the  54  hour  forecast.  In  it,  there  is  a  large  area  of  PVA  in  the  northern  Mediterranean 
between  Spain  and  Italy  that  is  much  stronger  for  the  control  run  than  for  Run3c.  The 
influence  of  this  PVA  can  be  seen  in  the  SLP  depiction.  Figure  26  shows  the  6  hour  prior 
SLP  and  the  SLP  for  the  concurrent  time  as  the  PVA.  Initially,  Run3c  pressure  contour 
was  only  slightly  above  the  control  run  pressure,  seen  as  a  slight  yellowing  over  the 
northern  Mediterranean.  As  the  PVA  develops  above  the  area,  a  low  pressure  system 
begins  to  initiate  for  the  control  run,  seen  as  the  addition  of  the  1000  mb  contour  in 
Figure  26.  Run3c  without  the  overlaying  PVA  maintains  the  same  SLP  level  seen  as  a 
deepening  of  yellow  (height  above  control)  in  the  same  area.  By  the  60  hour  forecast  the 
500  mb  PVA  is  lost  (figure  not  shown)  and  the  low  pressure  system  for  the  control  run 
loses  the  upper  level  support  and  therefore  fdls.  After  the  upper  PVA  diminished  Run3c 
SLP  heights  continue  to  be  greater  than  the  control  SLP  as  seen  in  Figure  27. 

b.  Case  4 

Figure  28  shows  the  comparison  of  the  thermal  advection  between  the 
control  run  and  Run3c  for  the  active  period  at  the  138  hour  forecast.  The  control  run  had 
a  more  developed  cyclonic  structure  at  700  mb  as  can  be  seen  by  the  deeper  low  pressure 
contours  over  Italy.  The  1000-500  mb  thickness  contours  appear  to  have  the  classic  self- 
development  structure  associated  with  cyclonic  development  at  the  surface.  Stronger 
thermal  advection  can  also  be  seen  for  the  control  run  by  the  deeper  shade  of  yellow  in 
eastern  Europe.  A  preexisting  surface  low  pressure  system  took  advantage  of  the  thermal 
advection  and  deepens  as  can  be  seen  by  the  extension  of  the  1000  mb  contour  within  the 
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core  of  the  cyclone  as  can  be  seen  in  Figures  29  and  30.  Run3c  maintains  the  pressure  of 
the  surrounding  1020  mb  pressure  surface  for  the  entire  period. 

3.  Central  Pressure  Low 

In  both  the  active  and  the  inactive  period,  there  was  a  notable  low  pressure 
development  near  the  center  of  the  domain  for  the  control  run.  There  was  little  support 
from  the  Petterssen  advection  terms  to  suggest  such  a  depression.  Although  the  previous 
paragraph  discussed  the  strengthening  of  the  central  low  pressure  system  by  thermal 
advection,  the  cyclone  was  already  well  developed  by  the  138  forecast  when  the  thermal 
advection  commenced. 

For  the  active  period,  the  control  run  low  pressure  system  began  to  develop  at  the 
102  hour  forecast  (figure  not  shown),  and  continues  to  develop  into  a  20  mb  depression 
by  the  end  of  the  forecast  seen  Figure  11.  The  location  of  the  depression  was  in  the 
center  of  the  inner  domain  over  Hungary.  A  similar  less  developed  system  for  Run3c  can 
be  seen  as  a  darkened  blue  area  over  western  Russia. 

In  the  inactive  period,  another  low  pressure  system  can  be  seen  for  the  control  run 
over  Italy.  Again  the  low  pressure  system  was  in  the  range  of  20  mb  and  was  located 
toward  the  center  of  the  model  domain.  For  this  system  the  development  began  at  about 
the  126  hour  forecast. 

A  check  was  made  to  determine  if  a  difference  in  terrain  height  was  contributing 
differences  in  the  central  domain  low  pressure  system.  During  the  COAMPS  Analysis 
program,  a  terrain  height  field  is  generated  that  determines  the  elevation  parameter  to  use 
for  each  grid  point,  based  on  the  latitude  and  longitude  positions.  As  discussed  earlier  in 
section  A-l,  the  lat/long  positions  do  not  match  up  perfectly  between  the  control  run  and 
each  of  the  test  runs.  The  reason  for  this  is  that  the  grid  point  positioning  is  first 
determined  based  on  the  coarse  grid  domain,  and  then  the  domain  spacing  is  reduced  by  a 
factor  of  three  for  the  higher  resolution  fine  grid  domain.  Since  the  comparison  grid 
points  were  first  based  on  the  coarse  grids,  which  were  different  (see  Figure  5),  the 
lat/long  position  of  the  fine  grid  points  were  slightly  different.  The  differences  in  the 
lat/long  positions  of  the  control  and  test  runs  were  small  (see  Table  1).  Still,  for 
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mountainous  regions,  the  slight  differences  in  position  induced  a  large  difference  in 
elevation  as  can  be  seen  in  Figure  31.  Figure  31  also  shows  that  there  is  not  a  “problem” 
with  the  terrain  height  in  the  central  area  of  the  domain  that  would  account  for  the  central 
pressure  low  that  develops  late  in  the  forecast  period.  The  average  difference  in  elevation 
over  the  entire  domain  was  less  than  a  meter. 

The  low  pressure  systems  discussed  here  could  not  be  supported  by  the  Petterssen 
Development  advection  terms  calculated  in  this  study.  Some  other  mechanism  drove  the 
cyclogenesis.  Run3c,  having  the  largest  coarse  nest  boundary  was  assumed  to  be  the 
model  with  the  closest  to  actual  observed  data,  since  it  would  have  the  least  adverse 
effects  from  the  lateral  boundary.  Therefore,  the  low  pressure  system  observed  in  the 
control  run  may  not  actually  exist.  The  fact  that  a  low  pressure  system  was  generated  for 
the  active  and  inactive  periods  in  the  same  general  location  at  nearly  the  same  forecast 
time  leads  to  questions  as  to  if  they  occurred  by  coincidence,  or  possibly  by  non-realistic 
forcing.  The  observed  control  run  central  domain  low  pressure  system  was  not 
investigated  further  in  this  study.  Additional  COAMPS  data  sets  would  be  required  to 
see  how  often  a  low  pressure  system  developed  in  the  latter  stages  of  an  extended 
forecast,  and  to  make  further  conjecture  as  to  potential  development  processes.  The  two 
periods  under  consideration  in  this  study  were  the  only  two  of  the  FNMOC  data  set 
carried  out  past  48  forecast  hours.  As  stated  previously,  a  special  request  was  made  of 
FNMOC  post  study  period,  to  extend  the  forecast  of  the  two  data  sets  under  consideration 
for  this  study. 

C.  CONCLUSIONS 

In  this  thesis,  an  experiment  was  devised  to  determine  the  significance  of  the 
separation  distance  between  the  inner  and  outer  nests  of  a  COAMPS  model  run.  Four 
COAMPS  runs  were  used  in  the  experiment.  The  first  data  set,  used  as  the  control,  was 
FNMOCs  European  Operational  Area  model.  The  other  three  data  sets  had  the  same  set¬ 
up  parameters  and  settings  as  the  control.  The  only  difference  was  the  size  of  the  outer 
boundary  area  and,  therefore,  the  outer-to-inner  nest  separation  distance.  There  were 
only  two  CO  AMPS  forecast  time  periods  used  in  this  study. 
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The  COAMPS  sea  level  pressure  fields  were  used  as  the  comparison  fields.  The 
choice  of  sea  level  pressure  was  determined  by  the  relevance  of  the  sea  level  pressure  to 
the  forecaster  and  more  specifically,  to  the  naval  operational  planner.  The  forecasting  of 
weather  conditions  at  the  surface  will  influence  the  operational  planner’s  decision  making 
process.  Forecasted  weather  conditions  are  a  vital  part  of  operational  planning. 

After  comparing  the  entire  domain  fields  between  the  control  and  the  test  fields,  it 
was  determined  that  observed  differences  in  the  models  were  slight.  The  control  run  used 
as  a  basis  field  contains  a  certain  degree  of  natural  variability.  When  the  RMS  difference 
between  the  control  run  and  the  test  runs  was  compared  to  the  natural  variability  of  the 
control  run,  it  was  found  that  the  differences  between  the  models  were  insignificant.  The 
sea  level  pressure  as  determined  by  the  four  model  runs  were  essentially  the  same  when 
the  entire  model  domain  was  compared  at  the  grid  point  level. 

When  a  direct  comparison  of  the  sea  level  pressure  fields  was  made  between  the 
two  model  runs  with  the  largest  difference  in  outer  domain  size,  isolated  areas  of 
significant  differences  were  observed.  Four  individual  cases  were  studied  further  to  see 
how  the  cyclonic  development  differed  due  solely  to  the  coarse  domain  size.  It  was 
found  that  the  COAMPS  boundary  domain  set  up  could  change  the  resultant  forecasted 
field  in  isolated  areas  that  would  change  the  predicted  cyclogenesis. 

This  study  only  looked  at  two  forecast  periods.  To  get  a  more  complete 
understanding  of  how  the  inner  and  outer  domain  setup  affects  the  resultant  COAMPS 
data  fields,  a  more  extensive  study  should  be  conducted  with  a  larger  set  of  data.  Also, 
only  the  sea  level  pressure  was  extensively  compared.  The  only  upper  level 
considerations  made  were  the  contributions  required  by  the  Petterssen  Development 
Equation.  Additional  study  could  be  made  on  comparisons  of  the  jet  or  the  500  mb 
contour  level. 

Another  practical  area  of  study  would  be  to  compare  the  four  data  sets  with 
observed  data  to  determine  which  of  the  four  performed  better,  and  if  a  smaller  outer 
boundary  is  “close  enough”  to  observed  weather  patterns. 
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The  other  area  of  study  mentioned  in  the  previous  section  would  be  to  fully 
explore  the  possible  bias  of  the  control  model  toward  the  development  of  a  low  pressure 
system  toward  the  end  of  the  forecast  period.  The  occurrence  of  the  low  pressure  system 
may  be  boundary  condition  related,  since  it  mainly  showed  up  in  the  control  run. 
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APPENDIX  A 


A.  SELECTED  FIGURES  PRESENTED  IN  THE  TEXT 

Active  Jet 


Figure  6  Jet  strength  determined  by  the  300  mb  isotachs  from  NOGAPS  Analysis  model 
data.  Units  are  in  meters  per  second.  Frames  are  given  in  12  hour  increments. 
Active  jet  shown  on  the  west  boundary  of  the  COAMPS  domain. 
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Figure  7  Jet  strength  determined  by  the  300  mb  isotachs  from  NOGAPS  Analysis  model 
data.  Units  are  in  meters  per  second.  Frames  are  given  in  12  hour  increments. 
Inactive  jet  period  shown  as  a  non  consistent  flow  on  the  western  boundary  of  the 
COAMPS  domain. 
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Figure  8  Sea  Level  Pressure  contours  for  the  active  period.  Forecast  hour  given  in  the 
lower  right  hand  corner.  24  hour  forecast  not  available  for  analysis. 
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Figure  9  Horizontal  depiction  for  the  active  period,  ‘000’  forecast  SLP  difference 
between  Run3c  and  the  control  run.  Only  slight  differences  of  less  than  5 
mb  are  present.  Control  run  contour  lines  are  given  in  black.  Blue  areas 
represent  areas  where  Run3c  contour  is  below  the  control  run. 
Red/Yellow  areas  represent  areas  where  Run3c  contour  is  above  the 
control  run. 
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Figure  10  Horizontal  depiction  for  the  active  period,  120  hour  forecast  SLP 
difference  between  Run3c  and  the  control  run.  Control  run  contour  lines 
are  given  in  black.  Blue  areas  represent  areas  where  Run3c  contour  is 
below  the  control  run.  Red/Yellow  areas  represent  areas  where  Run3c 
contour  is  above  the  control  run. 
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Figure  1 1  Horizontal  depiction  for  the  active  period,  144  hour  forecast  SLP 
difference  between  Run3c  and  the  control  run.  Control  run  contour  lines 
are  given  in  black.  Blue  areas  represent  areas  where  Run3c  contour  is 
below  the  control  run.  Red/Yellow  areas  represent  areas  where  Run3c 
contour  is  above  the  control  run. 
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Figure  12  Standard  deviation  of  the  control  run  sea  level  pressure  for  both  active  and 
inactive  periods. 
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Figure  13  Root  mean  square  difference  between  the  test  runs  and  the  control  run, 
active  period.  The  minimum  standard  deviation  of  the  control  run  is  also 
shown  for  reference.  RMS  of  Run3a  shown  in  blue.  RMS  of  run3b 
shown  in  green.  RMS  of  Run3c  shown  in  magenta. 
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Figure  14  STD  and  RMS  evaluation  of  domain  grid  points  less  than  990  mb,  active 
period.  Control  run  standard  deviation  is  given  in  red.  RMS  difference 
between  control  run  and  Run3a/b/c  are  given  in  blue,  green  and  magenta 
respectively. 
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Table  2  Percentage  of  subset  points  in  domain  after  correcting  for  grid  mismatch. 
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Control  Run 


Figure  15  Comparison  of  the  control  run  and  Run3c  Thermal  Advection  at  the  84  hr 
forecast,  active  period.  700  mb  contour  lines  shown  as  dashed  green  lines. 
500-1000  mb  thickness  line  shown  in  black.  Positive  thermal  advection 
shown  as  yellow/red  areas.  Negative  thermal  advection  shown  in  blue. 
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Figure  16  Horizontal  depiction  for  the  active  period,  78  and  84  hour  forecast  SLP 
difference  between  Run3c  and  the  control  run.  Control  run  contour  lines 
are  in  black.  Shaded  areas  represent  same  level  differences  as  before. 
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Control  Run 


Figure  17  Comparison  of  the  control  run  and  Run3c  Vorticity  Advection  at  the  96  hr 
forecast,  active  period.  500  mb  contour  lines  shown  in  black.  Red 
contours  indicate  positive  relative  vorticity.  Blue  contours  indicate 
negative  relative  vorticity.  The  contour  intervals  are  the  same  for  both  the 
positive  and  the  negative  vorticity.  Positive  relative  vorticity  advection 
shaded  in  yellow/red.  Negative  relative  vorticity  advection  shaded  in  blue. 
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Figure  18  Horizontal  depiction  for  the  active  period,  90  and  96  hour  forecast  SLP 
difference  between  Run3c  and  the  control  run.  Control  run  contour  lines 
are  in  black.  Shaded  areas  represent  same  level  differences  as  before. 
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Figure  19  Horizontal  depiction  for  the  active  period,  102  and  108  hour  forecast  SLP 
difference  between  Run3c  and  the  control  run.  Control  run  contour  lines 
are  in  black.  Shaded  areas  represent  same  level  differences  as  before. 
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Figure  20 


Sea  Level  Pressure  contours  for  the  inactive  period.  Forecast  hour  given  in  the 
lower  right  hand  comer.  24  hour  forecast  not  available  for  analysis. 
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Figure  21  Horizontal  depiction  for  the  inactive  period,  ‘000’  forecast  SLP  difference 
between  Run3c  and  the  control  run.  Only  slight  differences  of  less  than  5 
mb  are  present.  Control  run  contour  lines  are  given  in  black.  Blue  areas 
represent  areas  where  Run3c  contour  is  below  the  control  run. 
Red/Yellow  areas  represent  areas  where  Run3c  contour  is  above  the 
control  run. 
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Figure  22  Horizontal  depiction  for  the  inactive  period,  120  hour  forecast  SLP 
difference  between  Run3c  and  the  control  run.  Control  run  contour  lines 
are  given  in  black.  Blue  areas  represent  areas  where  Run3c  contour  is 
below  the  control  run.  Red/Yellow  areas  represent  areas  where  Run3c 
contour  is  above  the  control  run. 
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Figure  23  Root  mean  square  difference  between  the  test  runs  and  the  control  run, 
inactive  period.  The  minimum  standard  deviation  of  the  control  run  is 
also  shown  for  reference.  RMS  of  Run3a  shown  in  blue.  RMS  of  run3b 
shown  in  green.  RMS  of  Run3c  shown  in  magenta. 
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Figure  24  STD  and  RMS  evaluation  of  domain  grid  points  less  than  1010  mb, 
inactive  period.  Control  run  standard  deviation  is  given  in  red.  RMS 
difference  between  control  run  and  Run3a/b/c  are  given  in  blue,  green  and 
magenta  respectively. 
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Table  3  Percentage  of  subset  points  in  domain  after  correcting  for  grid  mismatch. 
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Control  Run 


Figure  25  Comparison  of  the  control  run  and  Run3c  Vorticity  Advection  at  the  54  hr 

forecast,  inactive  period.  500  mb  contour  lines  shown  in  black.  Red 
contours  indicate  positive  relative  vorticity.  Blue  contours  indicate 
negative  relative  vorticity.  Positive  relative  vorticity  advection  shown  as 
yellow/red  shaded  areas.  Negative  relative  vorticity  advection  shown  as 
shaded  blue  areas. 
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Figure  26  Horizontal  depiction  for  the  inactive  period,  48  and  54  hour  forecast  SLP 
difference  between  Run3c  and  the  control  run.  Control  run  contour  lines 
are  in  black.  Shading  represents  the  same  level  differences  as  before. 
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Figure  27  Horizontal  depiction  for  the  inactive  period,  60  and  66  hour  forecast  SLP 
difference  between  Run3c  and  the  control  run.  Control  run  contour  lines 
are  in  black.  Shading  represents  the  same  level  differences  as  before. 
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Control  Run 


Figure  28  Comparison  of  the  control  run  and  Run3c  Thermal  Advection  at  the  138  hr 
forecast,  inactive  period.  700  mb  contour  lines  shown  as  dashed  green 
lines.  500-1000  mb  thickness  line  shown  in  black.  Positive  thermal 
advection  shown  as  yellow/red  areas.  Negative  thermal  advection  shown 
in  blue. 
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Figure  29  Horizontal  depiction  for  the  inactive  period,  132  and  138  hour  forecast 
SLP  difference  between  Run3c  and  the  control  run.  Control  run  contour 
lines  are  in  black.  Shading  represents  the  same  level  differences  as  before. 
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Figure  30  Horizontal  depiction  for  the  inactive  period,  144  hour  forecast  SLP 
difference  between  Run3c  and  the  control  run.  Control  run  contour  lines 
are  in  black.  Shading  represents  the  same  level  differences  as  before. 
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Figure  31  Terrain  height  difference  between  control  run  and  Run3c  in  meters. 

Although  the  scale  is  only  given  for  up  to  20  meters  difference,  actual 
differences  in  isolated  areas  were  as  high  as  470  meter.  The  mean  height 
difference  was  0.38  meters  with  a  standard  deviation  of  12.8  meters. 
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APPENDIX  B 


A.  EQUATIONS  USED  IN  NOGAPS  MODEL  (FROM  HOGAN  AND 
ROSMOND,  1991) 

The  equations  for  the  NOGAPS  model  are  formulated  in  a  spherical  coordinates 
and  a  hybrid  vertical  coordinate.  The  horizontal  coordinates  are  the  longitude  A  and  the 
latitude#?.  The  vertical  coordinate  is  represented  by  //  and  has  a  range  from  0  to  1.0. 
From  the  principle  of  mass  continuity 


dn 

~dt 


u 


dp 

dp 


Ps 

drj-  |  V»(u  dp), 

Ptop 


(0.24) 


where  n  is  defined  as 


K  =  Ps-Ptop  (°-25) 

p  =  lmb  is  the  pressure  at  the  top  of  the  model;  ps  is  the  terrain  pressure;  and  dp  is  a 
function  of  A,  (p,  and  // .  The  vertical  motion  equation  is 
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where 


p{A,(p,p)-ptop 

s(';)-  MM)  ■ 


(0.27) 


The  thermodynamic  variable  used  by  the  forecast  model  is  the  virtual  potential 
temperature  0 ,  which  is  defined  by 


0  =  T 


^1  +  0.608 
P 


where  q  is  the  specific  humidity.  The  function  P  is  defined  as 


(0.28) 
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(0.29) 


P  = 


K 


V PoJ 


where  p0  is  1000  hPa  and  k  is  the  gas  constant  divided  by  the  heat  capacity,  k  =  R/cp. 
The  thermodynamic  equation  is  given  by 


d6_ 

dt 


rd6} 
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+  Qe> 
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(0.30) 


where  Qe  is  the  diabatic  heating.  The  quantity  Qe  is  the  sum  of  the  radiation,  latent  heat 

release,  and  the  horizontal  and  vertical  diffusion  processes.  The  advection  term  is 
defined  by 
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The  hydrostatic  equation  is  written  in  the  form 


d<f) 


(0.32) 


where  (j)  is  the  geopotential  and  P  is  given  by  (0.29).  The  above  form  is  chosen  to 
facilitate  total  energy  conservation  in  the  discrete  adiabatic  equation. 

For  the  NOGAPS  spectral  model,  the  choice  of  variables  used  to  describe  the 
atmospheric  motion  is  the  vorticity  £,  and  the  divergence  D.  The  £,  and  the  D  equations 
are  given  by 

A  =  (0.33) 

dt  cos  V  dA  dp 
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where  p ,  U,  V,  G,  H,  and  /  are  defined  as 


p  =  sin^, 


(0.35) 
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and 


U  =  u 


V  =  v 


cos  cp 
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(0.40) 


where  Qu  and  Qv  are  the  components  of  the  momentum  stress  and  Q  is  the  earth’s 
angular  rotation  rate. 

All  the  terms  in  the  equations  of  motion  are  written  in  spherical  harmonic  form 
using  the  transform 

i  l  2k 

X“(t)  =  —  1 1  X(/l ,ju,t)  x  P™ (jU^exp(-imA)dAdju,  (0.41) 

-1  o 

where  m  is  the  zonal  wavenumber  and  n  the  total  wavenumber. 

For  a  complete  description  of  the  NOGAPS  numerical  model  see  Hogan  and 
Rosemond,  1991. 
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APPENDIX  C 


A.  EQUATIONS  USED  IN  COAMPS  MODEL  (FROM  HODUR,  1997) 

The  COAMPS  atmospheric  model  is  comprised  of  the  nonhydrostatic, 
compressible  form  of  the  dynamics  and  parameterization  for  subgrid-scale  mixing, 
surface  fluxes,  explicit  moist  physics,  cumulus  convection  process  and  radiation.  The 
vertical  coordinate  <x  is  chosen  to  allow  for  flow  over  irregular  surfaces.  The  sigma 
coordinates  are  defined  as 


O’  =  ztop(' 


z-z. 


sfc 


"top 


"sfc 


(0.42) 


where  ztop  is  the  depth  of  the  model  domain  and  zsfc  is  the  terrain  height.  The  complete 
set  of  prognostic  equations  is  as  follows: 
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The  above  equations  were  determined  from  the  equation  of  state 


P  =  PRdTv 


the  definition  of  the  virtual  temperature 

rv  =  r(i.o+o.608<7v) 


and  the  Exner  function 
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For  the  above  equations,  p  is  the  pressure;  p  the  density;  Rd  the  gas  constant  for 
dry  air;  T  the  temperature;  qv,qc,  qr,  q,,  and  qs  the  mixing  ratios  of  water  vapor,  cloud 
droplets,  raindrops,  ice  crystals,  and  snowflakes,  respectively;  p00  a  constant  reference 
pressure;  Cp  the  specific  heat  at  constant  pressure  for  the  atmosphere;  u,  v,  and  w  the 

wind  components  in  the  x,  y,  and  z  directions,  respectively;  /  the  Coriolis  force;  g  the 
acceleration  due  to  gravity;  Sv,Sc,Sr,Sj,andSs  sources  and  sinks  of 

qv,qc,qr,  qt,  and  qs,  respectively;  Q0  sources  and  sinks  of  heat;  Vr  and  Vs  the  terminal 
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velocities  of  raindrops  and  snowflakes,  respectively;  0std  the  standard  atmospheric 
temperature;  qv  the  standard  mixing  ratio  corresponding  to  the  standard  atmospheric 
temperature;  and  D3  the  density  and  potential  temperature-weighted  three-dimensional 
divergence, 
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The  speed  of  sound  for  the  mean  state,  c  ,  is  defined  as 


c  = 


'  CAM''"2 

v  C.  y 


(0.62) 


where  Cv  is  the  specific  heat  at  constant  volume  for  the  atmosphere.  In  equations  (0.43)- 
(0.53)  and  (0.62),  the  prognostic  variables  n  and  6  have  been  decomposed  as 

n  =  n  +  7r' 

6  =  0  +  0'  (0.63) 


where  the  overbar  denotes  the  initial  mean  state,  which  is  a  function  of  z  only,  and  the 
prime  denotes  deviations  from  the  mean.  The  mean  state  follows  the  hydrostatic 
relationship 


dn 

da 


g 
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(0.64) 


For  a  complete  description  of  equations  used  in  the  CO  AMPS  model,  see  Hodur, 

1997. 
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APPENDIX  D 


HPC  DISTRIBUTED  CENTER  WORK-PLAN 

PROJECT:  COAMPS  OPERATIONAL  CONFIGURATION  TESTING 
I.  OBJECTIVE 


A  limited  area  model  (LAM)  is  formulated  so  as  to  obtain  temporal  changing  synoptic 
information  via  lateral  boundary  conditions  from  a  global  forecast  model.  The  temporal 
frequency  and  LAM  grids  over  which  the  energy  exchange  occurs  between  the  two  models  has  to 
be  fine-tuned  so  that  large-  and  short-wave  synoptic  features  from  the  global  model  propagate 
appropriately  to  the  LAM.  Currently,  FNMOC  runs  COAMPS  operationally  as  a  nested  LAM. 
The  lateral  boundary  conditions  to  the  COAMPS  outer  (coarse)  grid  are  provided  at  every  coarse 
grid  time-step  by  linearly  interpolating  global  NOGAPS  forecast  information,  which  are  input  for 
every  three  hour  COAMPS  forecast  interval.  Although  this  formulation  has  been  giving 
satisfactory  COAMPS  forecasts,  it  is  of  interest  to  determine  if  such  forecasts  can  be  improved 
by:  (i)  increasing  the  sampling  frequency  of  the  NOGAPS  fields  to  one  hour;  (ii)  changing  the 
relative  distance  between  COAMPS  inner  nest  and  its  coarse  grid  boundary  where  the  lateral 
boundary  conditions  are  introduced;  or  (iii)  running  the  COAMPS  forecasts  in  a  warm  start  mode 
to  reduce  phase  errors  for  synoptic  events. 


These  model  runs  will  be  conducted  to  determine  if  synoptic  meteorological  conditions 
(and  therefore  mesoscale  conditions)  are  better  posed  with  any  of  the  test  scenarios  outlined 
above  over  what  is  obtained  with  the  operational  setup.  NOGAPS  fields  run  specifically  to 
generate  lateral  boundaries  for  these  COAMPS  tests  will  be  used  as  ’’truth”  for  these  simulations 
as  it  is  assumed  that  one-way  COAMPS  forecasts  of  synoptic  conditions  can  be  no  more  accurate 
than  those  obtained  from  NOGAPS.  At  the  same  time,  a  fourth  (iv)  experiment  will  be  conducted 
whereby  the  outer  large  scale  COAMPS  nest  is  dropped  and  thus,  the  lateral  boundary  conditions 
from  NOGAPS  are  input  directly  to  the  current  inner  nest.  This  has  a  potential  cost  saving  and 
will  be  tested  to  determine  what,  if  any,  deterioration  in  the  forecast  is  observed  from  this  setup. 
A  fifth  COAMPS  configuration  will  test  the  impact  of  extending  the  current  operational  L30 
configuration  to  L40. 

It  is  proposed  that  six  parallel  numerical  simulations  be  completed  to  evaluate  each  of  the 
scenarios  mentioned  previously.  Verification  of  the  results  will  be  conducted  to  determine  the 
accuracy  of  both  synoptic  and  mesoscale  weather  patterns.  This  will  then  provide  a  first  step 
towards  optimizing  both  the  NOGAPS  sampling  frequency  and  COAMPS  operational 
configuration  with  the  goal  of  producing  the  most  accurate  forecasts  possible  within  given 
operational  constraints. 

II.  BACKGROUND 


From  long  experience  running  NOGAPS  operationally  at  FNMOC,  features  of  traveling 
waves  of  periods  as  small  as  an  hour  have  been  observed.  Such  short  wave  features  are  capable 
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of  impacting  weather  in  local  areas.  If  configured  properly,  COAMPS  would  be  expected  to 
resolve  such  weather  patterns.  However,  it  has  been  observed  that  COAMPS  forecasts 
occasionally  miss  entirely,  or  do  not  resolve  well  such  fast  moving  synoptic  short  wave  features. 
This  problem  is  common  to  all  LAM  forecasts  and  is  a  result  of  the  imperfect  treatment  of 
boundary  conditions  in  the  models. 

One  of  the  largest  sources  of  errors  in  using  boundary  information  for  the  LAM  is  the 
relatively  coarse  temporal  frequency  for  which  synoptic  information  is  transferred  to  the  LAM. 
Typically  at  FNMOC,  NOGAPS  forecast  boundary  information  is  transferred  to  COAMPS  at  a 
three  hour  forecast  interval.  However,  in  the  past,  NOGAPS  boundary  conditions  have  been 
provided  to  the  COAMPS  model  at  12  hour  intervals.  Through  a  set  of  sensitivity  experiments 
conducted  by  the  Naval  Research  Laboratory  (NRL),  significant  COAMPS  forecast 
improvements  were  observed  when  the  temporal  frequency  of  the  NOGAPS  boundary  conditions 
were  increased  to  three  hours.  However,  the  impact  of  increasing  these  boundary  conditions  to 
one  hour  has  never  been  explored. 

In  addition  to  the  temporal  frequency  issue,  there  has  been  interest  in  determining  if  the 
size  and  location  of  the  coarsest  grid  could  impact  the  accuracy  of  the  forecasts  produced  by  a 
LAM.  It  is  a  common  practice  for  LAM  models  to  employ  grid  structures  extending  as  far 
upstream  from  the  inner  grids  as  possible.  This  is  done  to  lessen  the  impact  of  the  imperfect 
boundary  condition  treatment  on  the  finer  scale  grids  interior  to  this  coarsest  grid.  However,  Cliff 
Mass  (University  of  Washington)  in  a  recent  discussion  with  FNMOC/NRL  suggested  that  the 
real-time  MM5  forecasts  conducted  at  the  University  of  Washington  show  marked  improvement 
in  the  synoptic  patterns  reproduced  by  MM5  when  its  coarsest  meteorological  MM5  grid  is 
located  as  close  to  the  interior  grids  as  possible. 

A  third  issue  to  be  explored  centers  around  the  use  of  warm  starting  COAMPS  forecasts 
every  watch  as  is  done  for  the  MM5  forecasts  at  AFWA.  The  potential  advantage  to  this  is  a 
reduction  in  synoptic  scale  phase  errors  due  to  the  initialization  of  the  large  scale  from  NOGAPS 
every  watch.  This  might  lead  to  better  forecasts  of  synoptic  scale  meteorological  conditions  but 
has  the  potential  disadvantage  of  a  loss  of  mesoscale  forecast  information. 

The  fourth  issue  to  be  examined  is  one  posed  by  NRL.  This  issue  is  a  cost  cutting 
measure  only  as  it  involves  the  elimination  of  the  coarsest  resolution  COAMPS  nest.  The 
elimination  of  the  coarsest  nest  would  lead  to  applying  the  NOGAPS  three  hour  boundary 
conditions  directly  to  a  finer  scale  nest.  The  downside  of  this  recommendation  would  be  the 
potential  degradation  in  the  forecast  due  to  the  mismatch  in  scale  between  the  NOGAPS 
boundary  conditions  being  input  and  the  scale  represented  by  an  interior  COAMPS  grid.  This 
becomes  more  pronounced  as  the  interior  COAMPS  grid  uses  higher  resolution  grids  as  it  will 
resolve  more  mesoscale  features  not  represented  by  the  large  scale  global  model. 

The  fifth  issue  to  be  explored  is  an  examination  of  the  impact  of  increasing  the  vertical 
resolution  within  the  COAMPS  forecasts  run  at  FNMOC.  This  is  a  continuation  of  previous 
work  in  which  the  number  of  vertical  grids  employed  by  the  model  was  extended  from  30  to  40 
levels.  This  has  been  shown  by  NRL  to  more  realistically  transport  wave  energy  aloft  in  the 
model.  However,  testing  by  FNMOC  for  cases  involving  summertime  conditions  over  CONUS 
indicated  that  there  was  little  impact  by  increasing  the  vertical  resolution  of  the  model.  However, 
the  lack  of  significant  synoptic  forcing  in  summertime  conditions  over  CONUS  might  have 
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skewed  the  results  for  these  cases.  Therefore,  these  sets  of  forecasts  will  be  used  to  examine  the 
impact  of  the  increased  resolution  for  wintertime  conditions. 


III.  PROJECT  PLAN: 

This  project  is  projected  to  take  approximately  6-9  months  to  complete  and 
involve  a  number  of  FNMOC  personnel.  Europe  has  been  chosen  as  the  CO  AMPS  area 
to  be  studied  due  to  its  importance  to  Navy  Fleet  activities,  the  diversity  and  frequency  of 
mesoscale  weather  events  observed,  and  the  large  number  of  observations  (which  can  be 
used  in  the  validation  of  the  results).  The  project  will  involve  two  phases  outlined  below. 
The  area  coverage  by  the  grid  used  in  this  project  will  closely  match  the  existing 
operational  European  coverage.  However,  the  test  areas  will  employ  a  higher  resolution 
54/1 8km  grid  structure  and  will  therefore,  not  match  exactly  the  geometries  used 
currently  in  operations.  The  choice  was  made  to  use  this  higher  resolution  domain  due  to 
FNMOC  plans  to  switch  all  COAMPS  area  meteorological  grids  to  this  structure  in  the 
near  future. 


Phase  I.  Daily  Model  Runs 


A  set  of  six  separate  meteorological  forecasts  will  be  run  to  evaluate  each  of  the 
scenarios  discussed  previously.  COAMPS  meteorological  forecasts  for  each  of  the  five 
cases  will  be  run  beginning  15  December  2003  and  continuing  through  15  March  2004. 
NOGAPS  boundary  conditions  will  be  provided  by  the  beta  NOGAPS  forecasts  running 
on  the  operational  AMS  1/2  machines.  Included  below  is  a  description  of  each  of  the 
daily  forecasts  to  be  run.  Each  of  the  forecasts  will  involve  four  forecasts  per  day.  Two 
of  the  runs  will  be  48  hours  long  (real  time  runs)  with  the  other  two  being  only  6  hours 
(off  time  runs).  All  CO  AMPS  runs  will  be  initialized  from  the  NOGAPS  fields  and  will 
be  run  in  an  unclassified  mode.  All  cases  will  also  involve  running  the  Multivariate 
Optimal  Interpolation  (MVOI)  model  to  combine  a  first  guess  with  available  observations 
from  the  operational  data  directories  to  provide  the  for  CO  AMPS  forecast  model  with  its 
initial  fields. 


Base  Case 

The  base  case  will  be  set  up  to  mimic  the  operational  run  currently  conducted  at 
FNMOC.  However,  as  mentioned  above,  the  horizontal  resolution  of  the  model  will  be 
increased  from  81/27/6km  to  54/1 8km.  Boundary  conditions  to  be  used  by  the  run  will 
be  obtained  from  the  three  hour  real  time  beta  NOGAPS  forecast.  Operational  OCARDS 
will  be  used  to  generate  all  fields  normally  produced  by  the  operational  run  with  these 
fields  being  stored  to  a  private  database.  This  base  case  will  be  run  in  a  "hot  start"  mode 
after  the  initial  start  up  of  the  model.  (The  term  "hot  start"  refers  to  the  creation  of  a  first 
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guess  field  used  by  the  forecast  coming  from  a  previous  six  hour  COAMPS  forecast  as 
opposed  to  obtaining  the  first  guess  from  NOGAPS  forecast  fields.) 


Experiment  2:  Increased  temporal  boundary  conditions 

The  set  up  of  this  mn  will  follow  the  base  case  mn  with  the  exception  of  use  of  additional 
boundary  condition  data  from  the  NOGAPS  model.  For  this  case,  hourly  boundary  conditions 
from  NOGAPS  will  be  used.  As  before,  the  meteorological  fields  from  this  forecast  will  be 
stored  to  a  separate  private  database  for  retention  and  comparison  at  a  later  date 

Experiment  3:  Coarse  Nest  Influence 

Four  experimental  mns  will  be  conducted  to  determine  how  the  lateral  separation 
between  the  inner  nest  and  its  coarse  grid  affect  the  COAMPS  forecast.  The  inner  nest  for  all 
mns  will  remain  unchanged  and  cover  a  constant  geographical  location.  The  coarse  grid  will  be 
increased  in  size  in  each  of  the  experiments,  allowing  for  a  doubling  of  the  lateral  separation  on 
the  upwind  side.  The  downwind  side  of  the  coarse  grid  (eastern  and  southern  boundaries)  will 
remain  the  same  in  each  experiment. 

Experiment  4:  Cold  Start 

This  case  will  involve  obtaining  a  first  guess  field  for  every  12  hour  off-time  forecast 
from  the  NOGAPS  forecast  fields.  Therefore,  each  real  time  48  hour  forecast  will  obtain  its  first 
guess  fields  from  a  six  hour  forecast  which  has  been  "cold  started”  from  NOGAPS  initial 
conditions.  All  other  elements  of  this  simulation  will  follow  the  base  case  scenario. 


Experiment  5:  Coarse  Grid  Elimination 


This  case  will  involve  the  use  of  an  18  grid  only.  Boundary  conditions  for  this  case  will 
again  be  developed  from  the  AMS  1/2  beta  NOGAPS  mns.  However,  the  boundary  conditions 
will  be  applied  directly  to  the  18km  grid.  The  application  of  these  boundary  conditions  to  this 
grid  will  be  identical  to  their  application  to  the  coarsest  grid  in  the  base  case  (e.g.  3  hour  updates 
linearly  interpolated  to  the  18km  grid  spacing).  All  other  elements  of  this  simulation  will  follow 
the  base  case  scenario. 


Experiment  6:  Increased  vertical  resolution 


The  daily  forecasts  for  this  case  will  involve  higher  resolution  in  the  vertical  than  the 
base  case.  Most  of  the  new  levels  will  be  included  in  the  middle  troposphere  with  only  a  single 
level  added  to  the  atmosphere  below  1500m.  All  other  elements  of  this  simulation  will  follow  the 
base  case  scenario. 
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Roger  S  will  take  the  lead  for  this  aspect  of  the  project  will  take  an  estimated  3 
months  to  complete. 


Phase  II.  Case  Study  Selection 

This  phase  will  take  two  months  to  complete  and  will  involve  the  selection  and  mnning 
of  case  study  days  to  be  modeled;  four  or  five  cases  will  be  analyzed.  The  choice  of  case  study 
periods  will  be  based  on  two  criteria: 

1.  Subjective  analysis  of  the  operational  CO  AMPS  Europe  gmeta  files  to 
identify  important  synoptic  and  mesoscale  weather  events  of  the  region  (e.g. 
Mistral,  intense  precipitation,  etc.); 

2.  Identification  of  significant  differences  between  runs  in  the  validation 
results  from  phase  I. 

These  two  lists  of  dates  will  be  developed  independently  and  compared  to  find  common 
dates  on  both  lists.  These  results  of  this  comparison  will  then  be  used  to  determine  the  particular 
case  study  days  to  examine  in  more  detail. 

Following  the  selection  of  the  case  study  days  to  be  examined,  the  COAMPS  model  and 
fields  from  each  of  the  Phase  I  experiments  will  be  used  to  generate  a  new  24  hour  forecast  for 
the  case  study  days  using  a  high  resolution  6km  inner  nest  not  mn  for  the  daily  cases.  This  inner 
nest  will  be  tailored  to  the  specific  area  of  interest  to  determine  if  higher  resolution  fields  over  the 
area  lead  to  a  better  forecast  of  the  event.  (Validation  will  again  utilize  standard  observations  and 
therefore,  grid  placement  will  also  be  dependent  on  the  availability  of  observations  in  an  area.) 

Jeff  L.  will  take  the  lead  in  determining  case  study  days  to  be  used  and  Roger  S.  will  be 
responsible  for  the  associated  model  mns. 


Phase  III.  Model  Validation  and  Verification 


There  are  two  levels  at  which  we  need  to  test  our  hypotheses:  Effect  of  increasing 
NOGAPS  sampling  frequency  and  the  effect  of  emplacing  the  inner  nest  closer  to  the  coarse-grid 
boundary  where  lateral  boundary  conditions  are  introduced,  on  the  (1)  propagation  of  the  short- 
and  long-wave  phenomenon  into  the  COAMPS  forecasts;  and  (2)  improvement  in  the  COAMPS 
forecasts. 


Measuring  the  improvement  in  a  mesoscale  model  forecast  involves  both  the 
subjective  and  objective  verification  of  sensible  weather  elements  that  are  routinely  measured 
or  estimated  (e.g.,  pressure,  temperature,  and  winds).  By  increasing  the  sampling  of  the 
lateral  boundary  conditions,  we  would  expect  to  see  an  overall  reduction  in  the  phase  error  of 
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synoptic  and  mesoscale  disturbances  entering  the  domain.  A  documented  subjective  analysis 
("man-in-the-loop")  of  the  selected  case  studies  will  be  carried  out  to  answer  questions  about 
phase  errors.  Of  equal  importance  and  to  complement  the  subjective  component,  an  objective 
verification  of  the  COAMPS  forecasts  against  station  observations  (or  satellite  surface  winds 
over  the  Mediterranean  Sea)  will  demonstrate  quantitatively  what,  if  any,  improvements  are 
observed  both  in  time  and  space.  Specifically,  objective  verification  will  involve: 

a)  Standard  verification  statistics  of  radiosonde  observations  for  Phase  I  results:  both 
height  and  time  series  errors.  Past  experience  running  "beta”  COAMPS  indicates  little 
departure  in  these  type  of  grid-averaged  statistics,  however,  this  information  will  be 
supplemented  with  other  parameters  and  methods  of  objective  verification  which  may  lead  to 
more  conclusive  results. 

b)  Surface  observation  verification  will  be  applied  to  both  Phase  I  and  II  results  and 
will  be  subsampled  to  zoom  in  on  the  area  of  interest  (Phase  II  results).  Parameters  to  be 
verified  include  winds,  temperature,  and  mslp. 

c)  Remotely  sensed  surface  wind  verification:  Standard  verification  scores  (e.g.,  bias, 
rmse)  will  be  calculated  against  SSM/I  and  Quikscat  scatterometer  maritime  winds.  Domain 
subsampling  may  reveal  systematic  differences  for  such  areas  as  Bay  of  Biscayne  or  the 
North  Sea  (increased  LBC  frequency  and  removal  of  nestl  experiments).  Event  verification 
software  (Nachamkin,  NRL)  may  also  reveal  systematic  model  experiment  differences. 
Contingency  table  statistics  will  also  be  computed  to  determine  any  improvement  in 
operationally  relevant  parameters  (Phase  I  and  II  results). 

d)  Raingauge  verification:  simple  differences  in  forecast-gauge  averaged  quantity.  If 
data  available,  It  is  anticipated  that  significant  experiment  differences  might  be  realized  with 
COAMPS  L40  forecasts  and  possibly  increased  LBC  frequency  in  the  form  of  better  wave 
depiction. 
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FNMOC  Distributed  Center 


Contributing  Sub-Project  Proposal 

Project  Title:  COAMPS  Lateral  Boundary  Conditions  Simulations 

Brief  Description:  Test  to  see  whether  incorporating  hourly  NOGAPS  boundary 
conditions  and/or  changes  to  the  operational  grid  structures  improve  COAMPS 
forecasts. 

Principal  Investigator:  Roger  Stocker 

Collaborator(s):  Andy  Hergert,  Torsten  Duffy,  Jeff  Lemer,  Ranjit  Passi,  Brad 
Harris  (NPS) 

Expected  start  and  stop  dates:  1  Dec  2003  /  30  June  2003 

Relate  this  project  to  DC  scientific  goals:  Improve  mesoscale  forecasting. 

What  operational  transitions  are  anticipated?  Should  the  results  be 
significantly  better  using  any  of  the  test  cases,  change  operational  COAMPS  to 
use  these  changes. 

Infrastructure  support  requirements: 

List  Unix  account  names  associated  with  this  sub-project:  ahergert, 
stockerr,  duffyt,  lemerj,  passir,  harrisb 
Typical  model  run  On-time  Off-time 

#  of  CPUs:  96  96 

Wall  clock  time:  2:00  :30 

#  of  runs  per  day:  16  16 

Access  classified  data?  No 

Request  high  priority  queue  (relative  to  other  DC  sub-projects)?  No 
Will  most  of  the  AMS3  resources  be  used  in  scheduled  (near  real  time)  model 
runs?  Yes 

If  yes, 

Recommended  Zulu  start  time:  +4,  with  forecasts  staggered  throughout 
the  watch 

Major  job  dependencies  (describe,  ex:  need  NOGAPS  boundary 
conditions):  NOGAPS  boundary  conditions  (hourly)  from  NOGAPS  beta 
run  on  AMS1.  This  will  require  code  changes  to  NOGAPS  (Duffy). 

File  system  volume  written  per  day,  indicate  MB/GB/TB 

/data  (indicate  'MODELS'  or  'NRL'):  30  GB  for  local  ISIS  files 
(MODELS) 

/scr:  120  GB  (MODELS) 

/tmp  (ops  only): 
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/  a /  ops/ etc/dynamic : 

ISIS  (indicate  alpha/beta/ops): 

File  system  volume  'running'  totals,  indicate  MB/GB/TB 

/data  (indicate  'MODELS'  or  'NRL'):  65  GB  (MODELS) 

/scr:  250  GB  (MODELS) 

/tmp  (ops  only): 

/  a /  ops/ etc/dynamic : 

ISIS  (indicate  alpha/beta/ops): 

Total  Project  Long  Term  Retention  (Silo) 

/data  (indicate  'MODELS'  or  'NRL',  MB/GB/TB):  5  TB  (MODELS) 
Retention  Period  (mos  /  yrs  /  indef):  6  mos 
(Note:  If  indefinite,  volume  should  be  per  year) 

Daily  transfers  OUT  of  AMS3  (internal  or  external): 

Volume  (MB/GB/TB): 

Number  of  files: 

Final  system  destination: 

Daily  transfers  IN  to  AMS3  (internal  or  external): 

Volume  (MB/GB/TB): 

Number  of  files: 

Originating  system: 

Additional  comments,  if  any: 
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APPENDIX  E 


A. 

DATAHD  FILES 

Runl 

Run3a 

Run3b 

Run3c 

16 

16 

16 

16 

total  number  of  pressure  levels 

30 

30 

30 

30 

total  number  of  full  sigma  levels 

2 

2 

2 

2 

map  proj.l:  mercator,2:  lambert,  3:  polar,  4:cartesian,  5:  spherical 

70 

70 

70 

70 

1st  latitue  where  map  projection  intersects  the  earth 

18 

18 

18 

18 

2nd  latitue  where  map  projection  intersects  the  earth 

12 

12 

12 

12 

long  aligned  north-south;  used  for  rotaing  lambert  or  polar  projec 

45 

48 

50 

51 

standard  latitude  used  to  set  up  coarse  domain 

12 

8 

4 

-2 

standard  longitude  used  to  set  up  coarse  domain 

54000 

54000 

54000 

54000 

coarse  domain  grid  spacing  (meter)  in  x  direction 

54000 

54000 

54000 

54000 

coarse  domain  grid  spacing  (meter)  in  y  direction 

2 

2 

2 

2 

maximum  number  of  domains  in  the  forecast 

7 

7 

7 

7 

real  data  LBC  identifier  for  nests.  6-PK,  7-Davis 

7 

7 

7 

7 

lateral  boundary  conditions  for  x  direction 

7 

7 

7 

7 

lateral  boundary  conditions  for  y  direction 

1 

1 

1 

1 

starting  i-coordinate  for  grid  stretching 

1 

1 

1 

1 

ending  i-coordinate  for  grid  stretching 

1 

1 

1 

1 

fractional  amount  of  stretching 

7 

7 

7 

7 

- UNKNOWN - 

121 

131 

141 

151 

number  of  grid  points  in  the  x  direction  for  coarse  domain 

88 

98 

108 

108 

number  of  grid  points  in  the  y  direction  for  coarse  domain 

1 

1 

1 

1 

mother  domain  i  points  that  defines  domain  1  staring  (1,1)  location 

1 

1 

1 

1 

mother  domain  j  points  that  defines  domain  1  staring  (1,1)  location 

61 

66 

71 

76 

coarse  domain  i  points  that  defines  coarse  domain  standard  lat 

44 

49 

54 

54 

coarse  domain  j  points  that  defines  coarse  domain  longitude 

1 

1 

1 

1 

- UNKNOWN - 

54000 

54000 

54000 

54000 

domain  1  i-grid  horizontal  resolution 

54000 

54000 

54000 

54000 

domain  1  j-grid  horizontal  resolution 

16 

15 

12 

10 

latitude  of  domain  1  staring  (1,1)  location 

0 

0 

0 

0 

longitude  of  domain  1  staring  (1,1)  location 

16 

16 

16 

18 

latitude  of  domain  1  lower  right  comer 

43 

43 

43 

43 

longitude  of  domain  1  lower  right  comer 

68 

73 

77 

78 

latitude  of  domain  1  ending  (m,n)  location 

360 

360 

360 

360 

longitude  of  domain  1  ending  (m,n)  location 

55 

55 

53 

47 

latitude  of  domain  1  upper  left  comer 

310 

297 

284 

279 

longitude  of  domain  1  upper  left  corner 

1 

1 

1 

1 

- UNKNOWN - 

1 

1 

1 

1 

- UNKNOWN - 

2 

2 

2 

2 

- UNKNOWN - 
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Runl 

Run3a 

Run3b 

Run3c 

295 

295 

295 

295 

number  of  grid  points  in  the  x  direction  for  domain  2 

196 

196 

196 

196 

number  of  grid  points  in  the  y  direction  for  domain  2 

12 

22 

32 

44 

mother  domain  i  points  that  defines  domain  2  staring  (1,1)  location 

12 

12 

12 

12 

mother  domain  j  points  that  defines  domain  2  staring  (1,1)  location 

148 

133 

118 

97 

domain  2  i  points  that  defines  centerdomain  latitude 

97 

112 

127 

127 

domain  2  j  points  that  defines  centerdomain  longitude 

1 

1 

1 

1 

- UNKNOWN - 

18000 

18000 

18000 

18000 

domain  2  i-grid  horizontal  resolution 

18000 

18000 

18000 

18000 

domain  2  j-grid  horizontal  resolution 

23 

23 

23 

24 

latitude  of  domain  2  staring  (1,1)  location 

0 

0 

0 

0 

longitude  of  domain  2  staring  (1,1)  location 

23 

23 

23 

24 

latitude  of  domain  2  lower  right  comer 

40 

40 

40 

40 

longitude  of  domain  2  lower  right  corner 

63 

63 

63 

64 

latitude  of  domain  2  ending  (m,n)  location 

360 

360 

360 

360 

longitude  of  domain  2  ending  (m,n)  location 

54 

54 

54 

55 

latitude  of  domain  2  upper  left  comer 

324 

324 

324 

323 

longitude  of  domain  2  upper  left  corner 

2 

2 

2 

2 

- UNKNOWN - 

7500 

7500 

7500 

7500 

vertical  grid  spacing  level  01 

5800 

5800 

5800 

5800 

vertical  grid  spacing  level  02 

4200 

4200 

4200 

4200 

vertical  grid  spacing  level  03 

2500 

2500 

2500 

2500 

vertical  grid  spacing  level  04 

1000 

1000 

1000 

1000 

vertical  grid  spacing  level  05 

1000 

1000 

1000 

1000 

vertical  grid  spacing  level  06 

750 

750 

750 

750 

vertical  grid  spacing  level  07 

750 

750 

750 

750 

vertical  grid  spacing  level  08 

750 

750 

750 

750 

vertical  grid  spacing  level  09 

750 

750 

750 

750 

vertical  grid  spacing  level  1 0 

750 

750 

750 

750 

vertical  grid  spacing  level  1 1 

750 

750 

750 

750 

vertical  grid  spacing  level  12 

1000 

1000 

1000 

1000 

vertical  grid  spacing  level  1 3 

1000 

1000 

1000 

1000 

vertical  grid  spacing  level  14 

1000 

1000 

1000 

1000 

vertical  grid  spacing  level  1 5 

1000 

1000 

1000 

1000 

vertical  grid  spacing  level  1 6 

800 

800 

800 

800 

vertical  grid  spacing  level  1 7 

800 

800 

800 

800 

vertical  grid  spacing  level  1 8 

800 

800 

800 

800 

vertical  grid  spacing  level  19 

600 

600 

600 

600 

vertical  grid  spacing  level  20 

400 

400 

400 

400 

vertical  grid  spacing  level  2 1 

300 

300 

300 

300 

vertical  grid  spacing  level  22 

200 

200 

200 

200 

vertical  grid  spacing  level  23 

140 

140 

140 

140 

vertical  grid  spacing  level  24 

90 

90 

90 

90 

vertical  grid  spacing  level  25 

60 

60 

60 

60 

vertical  grid  spacing  level  26 
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Runl 

Run3a 

Run3b 

Run3c 

40 

40 

40 

40 

vertical  grid  spacing  level  27 

30 

30 

30 

30 

vertical  grid  spacing  level  28 

20 

20 

20 

20 

vertical  grid  spacing  level  29 

20 

20 

20 

20 

vertical  grid  spacing  level  30 

31050 

31050 

31050 

31050 

sigmma  level  (meter)  01 

24400 

24400 

24400 

24400 

sigmma  level  (meter)  02 

19400 

19400 

19400 

19400 

sigmma  level  (meter)  03 

16050 

16050 

16050 

16050 

sigmma  level  (meter)  04 

14300 

14300 

14300 

14300 

sigmma  level  (meter)  05 

13300 

13300 

13300 

13300 

sigmma  level  (meter)  06 

12425 

12425 

12425 

12425 

sigmma  level  (meter)  07 

11675 

11675 

11675 

11675 

sigmma  level  (meter)  08 

10925 

10925 

10925 

10925 

sigmma  level  (meter)  09 

10175 

10175 

10175 

10175 

sigmma  level  (meter)  1 0 

9425 

9425 

9425 

9425 

sigmma  level  (meter)  1 1 

8675 

8675 

8675 

8675 

sigmma  level  (meter)  12 

7800 

7800 

7800 

7800 

sigmma  level  (meter)  13 

6800 

6800 

6800 

6800 

sigmma  level  (meter)  14 

5800 

5800 

5800 

5800 

sigmma  level  (meter)  15 

4800 

4800 

4800 

4800 

sigmma  level  (meter)  16 

3900 

3900 

3900 

3900 

sigmma  level  (meter)  17 

3100 

3100 

3100 

3100 

sigmma  level  (meter)  1 8 

2300 

2300 

2300 

2300 

sigmma  level  (meter)  19 

1600 

1600 

1600 

1600 

sigmma  level  (meter)  20 

1100 

1100 

1100 

1100 

sigmma  level  (meter)  21 

750 

750 

750 

750 

sigmma  level  (meter)  22 

500 

500 

500 

500 

sigmma  level  (meter)  23 

330 

330 

330 

330 

sigmma  level  (meter)  24 

215 

215 

215 

215 

sigmma  level  (meter)  25 

140 

140 

140 

140 

sigmma  level  (meter)  26 

90 

90 

90 

90 

sigmma  level  (meter)  27 

55 

55 

55 

55 

sigmma  level  (meter)  28 

30 

30 

30 

30 

sigmma  level  (meter)  29 

10 

10 

10 

10 

sigmma  level  (meter)  30 

10 

10 

10 

10 

pressure  level  (mb)  01 

20 

20 

20 

20 

pressure  level  (mb)  02 

30 

30 

30 

30 

pressure  level  (mb)  03 

50 

50 

50 

50 

pressure  level  (mb)  04 

70 

70 

70 

70 

pressure  level  (mb)  05 

100 

100 

100 

100 

pressure  level  (mb)  06 

150 

150 

150 

150 

pressure  level  (mb)  07 

200 

200 

200 

200 

pressure  level  (mb)  08 

250 

250 

250 

250 

pressure  level  (mb)  09 

300 

300 

300 

300 

pressure  level  (mb)  10 
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Runl 

Run3a 

Run3b 

Run3c 

400 

400 

400 

400 

500 

500 

500 

500 

700 

700 

700 

700 

850 

850 

850 

850 

925 

925 

925 

925 

1000 

1000 

1000 

1000 

31055 

31055 

31055 

31055 

26482 

26482 

26482 

26482 

23849 

23849 

23849 

23849 

20576 

20576 

20576 

20576 

18442 

18442 

18442 

18442 

16180 

16180 

16180 

16180 

13609 

13609 

13609 

13609 

11784 

11784 

11784 

11784 

10363 

10363 

10363 

10363 

9164 

9164 

9164 

9164 

7186 

7186 

7186 

7186 

5575 

5575 

5575 

5575 

3012 

3012 

3012 

3012 

1457 

1457 

1457 

1457 

762 

762 

762 

762 

111 

111 

111 

111 

pressure  level  (mb)  1 1 
pressure  level  (mb)  12 
pressure  level  (mb)  13 
pressure  level  (mb)  14 
pressure  level  (mb)  15 
pressure  level  (mb)  1 6 
standard  atmospheric  height,  level  01 
standard  atmospheric  height,  level  02 
standard  atmospheric  height,  level  03 
standard  atmospheric  height,  level  04 
standard  atmospheric  height,  level  05 
standard  atmospheric  height,  level  06 
standard  atmospheric  height,  level  07 
standard  atmospheric  height,  level  08 
standard  atmospheric  height,  level  09 
standard  atmospheric  height,  level  10 
standard  atmospheric  height,  level  1 1 
standard  atmospheric  height,  level  12 
standard  atmospheric  height,  level  13 
standard  atmospheric  height,  level  14 
standard  atmospheric  height,  level  15 
standard  atmospheric  height,  level  16 
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